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SUMMARY 

\ 

\ 

Like  the  free  propeller  in  axial  flight ,~a  single-rotation  ducted 
fan  of  highest  induced  efficiency  is  characterized  by  an  ultimate  wake 
vortex  system  shed  from  the  blade  trailing  edges  whose  apparent  motion 
is  that  of  rigid  helical  surfaces.  In  addition  and  concentric  with  this 
inner  sheet  system  there  is  a  cylindrical  surface  of  helical  vortex 
filaments  shed  from  the  duct  trailing  edge.  For  zero  hub  diameter  and 
neglecting  compressibixity,  viscosity,  and  tip  clearance,  a  consistent 
mathematical  model  of  the  constant -diameter  vortex  wake  is  developed 
and  the  compatibility  relationships  to  be  satisfied  are  presented. 

Using  the  Biot-Savart  equation,  the  vortex  strength  distribution  in 
the  ultimate  wake  is  determined  and  then  related  to  the  blade  bound 
vortex  strength  distribution.  In  addition,  expressions  are  developed 
for  the  thrust,  power,  and  induced  efficiency  which  depend  on  numerical 
integrations  of  velocity  profiles  in  the  ultimate  wake.  It  is  shown 
that  the  wake  vorticity  and  the  velocity  distribution  in  the  wake,  for 
all  loadings  from  the  lightly  loaded  limit  to  the  heavily  loaded  static 
thrust  condition  may  be  extracted  from  the  lightly  loaded  result  through 
the  use  of  a  simple  algebraic  scale  factor.  The  performance  parameters, 
thrust,  power  and  blade  bound  vortex  strength,  are  thus  expressed  in 
terms  of  a  lightly  loaded  solution  for  a  given  freestream  velocity, 
blade  tip  speed,  blade  number  and  the  loading  parameter.  The  lightly 
loaded  case  is  compared  to  existing  theoretical  and  experimental 
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results.  Some  sample  results  for  heavily  loaded  performance  parameters 
and  design  parameters  are  presented  and  compared  to  the  exact  results 
for  a  heavily  loaded  ducted  fan  with  an  infinite  number  of  blades. 
Results  for  the  heavily  loaded  system  over  a  broad  range  of  conditions 
are  presented.  ' 
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CHAPTER  I 

INTRODUCTION 

In  recent  years  there  has  been  an  increasing  interest  in  low 

speed  thrust  and  high-lift  devices  for  application  to  aircraft  with 

short  or  vertical  take-off  and  landing  characteristics.  One  of  the 

devices  under  consideration  for  this  role  is  the  ducted  fan  or  shrouded 

propeller.  Analytical  and  experimental  work  in  this  field  have  been 

summarized  in  an  excellent  review  of  the  state  of  the  art  through  1963 

by  Sacks  and  Burnell1.  As  pointed  out  in  the  review,  there  have  existed 

for  some  time  adequate  means  of  sizing  and  designing  a  duct  for  a  given 

choice  of  fan  blade  design  but  there  has  been  no  method  available  for 

optimizing  the  blade  design  for  the  finite  bladed  system  which  operates 

at  any  but  the  lightest  loadings.  For  the  lightly  loaded  case  the 

2  3  4  b 

problem  has  been  solved  by  Tachmindji  ,  Theodorsen  ,  and  Gray  *  .  In 

6  n  8 

more  recent  years  Morgan  ,  Ordway,  Sluyter,  and  Sonnerup',  Chaplin  , 

and  others  have  published  papers  in  the  field  of  ducted  fans,  still 

with  primary  emphasis  on  the  shroud  design.  It  is  the  purpose  of  this 

research  to  provide  the  information  needed  to  size  and  design  the 

optimum  heavily  loaded  ducted  fan. 

The  casis  for  this  work  stems  from  the  classical  analysis  of 

0 

the  free  propeller.  It  has  been  shown  by  Betz  that  an  isolated 
propeller  having  the  highest  possible  induced  efficiency,  that  is, 
an  optimum  free  propeller,  will  generate  an  ultimate  wake  vortex 
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system  which  moves  through  the  fluid  medium  as  if  the  vortex  sheets 

of  the  wake  formed  a  rigid  helical  structure  of  constant  pitch. 

Application  of  this  constraint  to  the  motion  of  the  sheets  provides 

a  straightforward  means  of  determining  the  radial  distribution  of  shed 

vortex  sheet  strength  distribution.  This  determination  has  been  carried 

out  by  Goldstein^  and  Theodor sen^. 

k 

It  has  been  shown  by  Gray  that,  for  the  ducted  fan  having  the 
highest  possible  induced  efficiency,  the  same  arguments  and 
considerations  are  valid  concerning  the  geometry  and  motion  of  those 
vortex  sheets  which  are  shed  from  the  blade  trailing  edges.  These 
arguments  are  reproduced  in  Appendix  I  for  completeness. 

In  addition  to  this  constraint,  the  system  must  satisfy  the 
Kutta  condition  at  the  duct  trailing  edge;  that  is,  the  flow  at  the 
trailing  edge  must  be  tangent  to  the  duct  mean  camber  surface.  As  a 
consequence  of  this  condition  a  sheet  of  vorticity  must  be  shed  from 
the  trailing  edge  forming  a  boundary  sheet  enclosing  the  screw-like 
sheets  shed  from  the  blade  trailing  edges.  It  is  the  determination  of 
the  geometry,  me  cion,  and  strength  distribution  of  this  boundary  sheet 
and  its  mutual  interaction  with  the  inner  sheets  which  constitute  the 
development  of  the  ultimate  wake  vortex  model  for  the  ducted  fan. 

When  the  geometry  and  motion  of  the  vortex  sheets  of  the  ultimate 
wake  have  been  defined,  a  straightforward  application  of  the  Biot-Savart 
equation  determines  the  strength  distribution  of  the  vortex  sheets. 

The  blade  bound  circulation  or  vortex  strength  distribution  is  then 
obtained  by  integration  of  the  inner  helical  sheet  strength.  This 
bound  circulation  is  the  basic  information  from  which  the  fan  blades 
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can  be  designed  using  any  of  the  methods  available  in  the  literature 
(e.g.  Theodorsen'1'1).  In  addition,  the  vortex  sheet  strength  distribu¬ 
tions  enable  a  detailed  calculation  of  the  flow  in  the  wake  from  which 


II 


the  thrust  and  induced  power  may  be  calculated. 

In  the  development  and  analysis  which  follows  the  effects  of 
hub  diameter,  blade  tip  clearance,  compressibility  and  viscosity  are 
neglected.  The  results  will  therefore  give  an  upper  limit  on  the  per¬ 
formance  of  single-rotation  fans  without  stator  vanes  for  comparison 
with  other  design  methods. 
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CHAPTER  II 

DEVELOPMENT  OF  THE  WAKE  VORTEX  MODEL 

Earlier  analyses  by  Gray  (see  Appendix  I)  have  shown  that  the 
optimum  condition  is  obtained  for  the  ducted  fan  when  the  helical 
vortex  sheets  shed  from  the  blades  appear  to  move  as  rigid  screw 
surfaces  in  the  ultimate  wake.  The  principle  difference  between  the 
arguments  used  by  Gray  and  those  for  the  free  propeller  is  that  the 
induced  velocity  in  the  wake,  at  the  surface  of  a  helical  vortex  sheet, 
need  not  be  normal  to  the  sheet  surface.  Beyond  the  need  for  satisfying 
the  Kutta  condition  at  the  duct  trailing  edge  no  other  information  is 
obtained  and  additional  relationships  must  be  developed  in  order  to 
assure  a  compatible  vortex  model. 

Further  information  may  be  obtained  through  consideration  of 
the  flow  associated  with  the  vortex  systems.  A  blade  trailing  vortex 
sheet  has  large  radial  velocity  components  associated  with  it  and,  with 
respect  to  a  coordinate  system  fixed  in  the  duct,  this  flow  is  periodic. 
Consequently,  if  the  duct  mean  camber  surface  is  to  correspond  to  a 
streamline,  a  distribution  of  vorticity  must  be  placed  along  this  surface 
to  cancel  all  normal  velocity  components.  Part  of  this  distribution 
must  be  periodic  in  nature  and  must  also  be  a  function  of  the  number  of 
blades  in  the  fan  and  must  rotate  with  the  blades.  A  portion  of  this 
periodic  distribution  is  considered  to  be  a  continuation  oi  each  blade 
bound  vortex  which  passes  directly  from  the  blade  tip  onto  the  duct 
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contour  where  it  spreads  out  over  the  duct  mean  camber  surface.  Then, 
in  accordance  with  the  Helmholtz  theorems,  it  moves  aft  along  the  mean 
camber  surface  and  is  shed  at  the  duct  trailing  edge  as  free  vortex 
filaments  of  varying  geometry  and  density  wrapped  on  a  cylindrical 
surface.  While  these  filaments  are  within  the  duct  contour  it  is 
assumed  that  they  make  all  of  the  necessary  adjustments  in  phase 
relationship  and  density  with  the  blade  trailing  vortex  system  so 
that  the  motion  of  this  wake  boundary  system  is  along  the  tangent  to 
the  mean  camber  surface  at  the  trailing  edge  of  the  duct  as  shown  in 
Pig.  1.  In  addition  to  the  blade  bound  vortex  continuation,  a  portion 
of  the  duct  trailing  vortex  sheet  is  considered  to  be  shed  from  bound 
vortices  of  varying  strength  rotating  in  the  duct  with  the  blades. 

It  is  also  assumed  that  the  duct  is  designed  so  that  there  is  no 
contraction  or  expansion  of  the  wake  downstream  of  the  fan.  This  is 

Q 

a  feasible  design  problem0  and  guarantees  that  the  duct  is  compatible 
with  the  wake  geometry  of  the  analysis. 

Continuing  with  the  examination  of  the  velocities  associated  with 
the  wake  vortex  sheets,  consider  a  helical  coordinate  system  r,  §,  £ 
defined  in  terms  of  the  cylindrical  coordinates  r,  Y ,  z  as  shown  in 
Fig.  2.  At  a  given  instant  an  inner  helical  vortex  sheet  in  the 
ultimate  wake  coincides  with  the  £  =  0  surface.  Then 


r  =  r 


0  s  r  j£  os 


%  =  tY  cos  cp  +  z  cos  cp  ,  ■  ®  s  5  s 


-  (V..  +  w)  COS  CD 

uu 

2b  (fi/ 2ft) 


(V  +  w)  cos 

■  CO 

2b  (n/2n) 


£  =  -  rf  sin  cd  +  z  cos  cp 


*  C  * 
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where  £  is  measured  along  the  helical  vortex  filaments  and  £  is  measured 
along  a  helical  line  normal  to  the  vortex  filaments.  It  should  be  noted 
at  the  outset  that  w,  the  apparent  velocity  parameter,  is  not  an  actual 
disturbance  velocity  associated  with  the  vortex  sheets  of  the  wake. 
Rather,  it  is  the  speed  along  the  wake  axis  with  which  the  inner  helical 
sheets  appear  to  be  moving  relative  to  a  coordinate  system  fixed  in  space. 
Similarly  (V^  +  w)  is  the  axial  speed  with  which  the  inner  sheets  appear 
to  move  relative  to  the  ducted  fan. 

If  the  apparent  rigid  structure  of  the  inner  helical  sheets  in 
the  ultimate  wake  is  to  be  maintained,  an  observer  fixed  to  any  point 
on  a  sheet  must  see  the  same  distribution  of  vorticity  in  the  wake 
regardless  of  his  angular  position  on  the  sheet.  He  must  further  see 
the  same  distribution  above  as  he  sees  below.  These  constraints  on 
the  distribution  of  vorticity  preclude  the  possibility  of  radial  dis¬ 
tortion  of  the  inner  sheets  and  amount  to  a  requirement  of  helical 
symmetry  of  the  wake  vortex  system.  Subject  to  this  helical  symmetry 
the  disturbance  velocity  vector  will  be  constant  along  the  helical 
lines  r  =  constant  and  £  =  constant  ooth  inside  and  outside  the  wake. 

Now,  consider  a  line  integral  of  the  velocity  along  the  path 
ABCDA  within  the  ultimate  wake  eis  shown  in  Fig.  3*  The  velocity  diagram 
with  respect  to  the  rotating  blades  is  shown  in  Fig.  4.  Along  BC  u^ 
is  constant  since  r  and  £  are  constant.  AB  and  CD  are  radial  lines 
at  two  values  of  ¥  on  the  sheet  surface  and  DA  coincides  with  the  wake 
axis.  Since  no  vorticity  is  enclosed  by  the  path,  the  line  integral 


1 


is  zero  so  that 


Velocity  Diagram  at  an  Inner  Helical  Vortex  Filament  of  the  Ultimate 
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pB  pC  pD  pA 

j  ur  to  +  J  \  +  J  ur  dr  +  J  d§  =  0 


B 


or 


^  (5C  "  ^  +  \  (ZA  -  ZD}  =  0 
=o 

where  u^  is  the  value  of  u^  on  the  axis  (r  =  0). 
Bnploying  the  helical  coordinate  relationships 


2tt 


V  +  w 
00 


-  u5  (r  T  cos  'P  +  b(57 aTT  3111 


in  p)  +  u. 


V  +  w 
00 

SQ  b(Q/2rr) 


=  0 


Then 


f  <V„  +  »)/®  , 

"s =  v  k _  t.  ♦ « — ) 


O  ^  COS  9  +  C° 


OR 


sin 


<P 


but  tan  9r  =  (Vm  +  w)/flR  and  for  the  coordinate  £ 
r  tan  9  =  R  tan  9p  so  that 


f _ tan  cp _ 1 

§  ~  u§  "Icos  9  +  tan  9  sin  9/ 


or 


u.  =  u-  sin  9 
s  So 

Consider  next  the  line  integral  along  the  path  EFGHE  as  shown 
in  Fig.  3-  The  portion  IS  is  at  r  =  R  .  From  the  requirement  of 
helical  symmetry,  the  filaments  of  the  boundary  vortex  sheet  must 
cross  the  line  of  intersection  with  an  inner  helical  sheet  at  every 
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point  on  the  intersection  line  at  the  same  filament  pitch  angle 
regardless  of  angular  position.  In  particular,  if  the  pitch  angle 
of  the  boundary  sheet  filaments  is  identical  to  that  of  the  inner 
sheet  as  the  boundary  is  reached  (r  =  R)  then  no  vorticity  is  enclosed 
by  the  path  so  that  the  line  integral  is  zero.  (The  equal  pitch 
constraint  is  a  lightly  loaded  condition  vhich  will  be  relaxed  later). 
Thus 


u?  (SG  -  5p)  +  u  sin  cpR  (SE  -  SH)  =  0 
or 


U_  =  u  sin  cp 
s  so 

outside  the  wake  as  well  as  inside. 

Further,  the  path  ABCDA  may  be  shifted  axially  so  that  the  helical 
portion  of  the  path  lies  along  an  arbitrary  helical  line  at  the  pitch 
angle  co  anywhere  inside  the  wake.  The  restrictions  of  symmetry  still 
require  u_  to  be  constant  along  this  path  so  that  u_  =  U-  sin  cp 
throughout  the  wake. 

These  results  are  subject  to  the  additional  constraint  that  the 
line  integral  of  velocity  taken  along  a  path  enclosing  the  wake  must 
be  zero.  The  line  integral  may  then  be  taken  along  the  path  ABCA  in 
Fig.  5  where  u^  along  ABC  is  constant  and  known  so  that 

u  dz  +  2rrR  u_  tan  co-  =  0 
<1  Z  |  R 

C  ° 

The  integral  is  then  taken  along  anpFA.  Subject  to  helical  symmetry 


l4 


the  induced  velocity  is  zero  along  the  radial  lines  CD  and  EA.  If  DE 

is  allowed  to  approach  infinity  then  u  along  DE  is  zero.  Since  no 

z 

vorticity  is  enclosed  by  the  path 

uz  dz  =  0 
C 


Then,  refering  to  the  previous  integration  and  subject  to  the  condition 

of  identical  filament  pitch  angle  along  the  intersection  line,  this 

result  requires  that  u_  be  zero.  Thus,  the  induced  velocity  must  be 

normal  to  the  inner  helical  sheets.  These  two  coexistent  conditions 

of  equal  filament  pitch  and  normal  induced  velocity  characterize  the 

4 

lightly  loaded  conditions  as  discussed  by  Gray  . 

The  more  general  condition  of  load  requires  only  that  the 
filaments  of  the  boundary  sheet  all  cross  the  intersection  line  at 
the  same  pitch  angle  so  that  a  relative  motion  between  the  inner  and 


outer  systems  may  exist.  Thus,  the  path  EFGHE  in  Fig.  3  will  enclose 
vorticity  and  no  information  will  be  obtained  from  the  integral  for 
a  heavily  loaded  ducted  fan.  However,  the  requirements  on  u^  still 
apply  inside  the  wake,  and  for  the  apparent  rigid  motion  of  the  inner 
helical  sheets  u^  must  be  proportional  to  cos  cp.  That  is,  from  the 
previous  results  and  the  velocity  diagram  of  Fig.  4 


u_  =  u,  sm  cp 
*  ^o 


=  w  cos  cp 


A  geometry  and  motion  must  now  be  established  for  the  boundary 
sheet  which  will  maintain  helical  symmetry  of  vortex  strength  distri¬ 
bution  and  which  will  permit  the  boundary  sheet  to  have  an  axial  motion 
relative  to  the  inner  sheets.  Further,  the  disturbance  velocities  must 
be  zero  outside  the  wake  in  order  that  the  flow  be  irrotational  there. 
This  may  be  shown  by  considering  the  line  integral  about  ACDEA  as  shown 

in  Fig.  5.  EA  and  CD  are  radial  lines  along  which  u  must  be  zero  by 

r 

helical  symmetry.  ED  lies  at  r  =  «  so  that  u  is  zero  along  this  path. 

z 

Since  the  path  encloses  no  vorticity  the  integral  of  u  along  AC  must 

z 

be  zero.  For  this  to  occur  either  u  is  identically  zero  or  u  is 

Z  L 

part  positive  and  part  negative.  From  considerations  of  continuity, 

for  the  latter  condition  to  occur  the  flow  must  form  closed  streamlines 

within  the  area  defined  by  the  path  of  integration.  However,  integration 

of  velocity  along  these  closed  streamlines  would  yield  finite  values  of 

circulation.  Thus,  a  contradiction  is  arrived  at  in  terms  of  the 

irrotationality  of  the  flow  so  that  ug  must  be  zero  along  AC.  Another 

line  integral  is  taken  along  ABCDEA.  For  the  helical  portion  of  the 

path  the  velocity  u^  must  be  constant  due  to  helical  symmetry.  Since 

the  velocities  along  the  other  branches  of  the  path  are  zero  and  no 

net  vorticity  is  enclosed,  u^  must  be  identically  zero  along  ABC. 

The  path  ABC  may  be  shifted  along  the  z-axis  to  any  axial  position 

with  no  change  in  the  results.  The  path  is  closed  by  joining  the 

helical  part  to  the  radial  parts  along  the  z-path  AC.  Thus  with 

u- ,  u  and  u  all  identically  zero  on  the  outer  surface  it  may  be 
5  z  r 

concluded  that  all  induced  velocities  arc  zero  on  the  outer  surface 
of  the  wake  boundary.  This  results  may  be  extended  to  include  the 
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entire  region  outside  the  wake  by  allowing  the  paths  of  integration 
to  expand  radially. 

The  strength  of  a  vortex  sheet  is  equal  to  the  discontinuity 
in  the  velocity  components  as  the  sheet  is  crossed  and  the  motion  of 
the  sheet  along  the  line  of  discontinuity  is  equal  to  the  mean  value 
of  the  velocity  at  R  and  R  .  At  r  =  R  ,  the  induced  velocity  at  the 
inner  surface  of  the  boundary  sheet  is  given  by  the  components 

u^  =  u^  sin  cpR 


while  the  induced  velocity  at  R+  is  zero.  Thus  at  the  line  of 
intersection  with  an  inner  helical  sheet  the  boundary  vortex  sheet 
strength  is  given  bjr 

y(Cb)  =  sin2  cpR  +  w2  cos2  cpR)2  (l) 


The  velocity  of  the  filaments  of  the  boundary  sheet  as  they  cross  the 
intersection  line  must  be  in  the  direction  of  the  induced  velocity 
at  R  and  normal  to  the  filament  direction;  that  is 


i/2  .2  ^  2 

ur  =  Sln  cp_  +  w 

5o  K 


2 

cos 


(2) 


The  filaments  of  the  boundary  vortex  sheet  must  then  all  cross 
the  lines  of  intersection  between  the  inner  screw  surface  and  the 
fivlind-rinal  hmindarv  at  a  rnnatanr.  an^le  m  not  equal  to  -■  This 

”  v  -  <-*  TR 

angle  may  be  detcrminevl  from  the  flight  speed  (V  ),  the  blade  rotational 
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speed  (OR)  and  the  total  disturbance  velocity  of  equation  (l).  The 
two  vortex  systems  are  related  through  equations  (l)  and  (2)  but  only 
along  their  lines  of  intersection.  On  the  wake  boundary  between  these 
helical  intersection  lines,  the  filament  density  or  sheet  strength  and 
the  filament  pitch  angle  will  vary  with  the  helical  coordinate  £. 

The  boundary  sheet  serves  several  basic  purposes.  First,  it 
must  cancel  the  radial  velocity  field  at  the  boundary  associated  with 
the  inner  helical  sheets.  Second,  it  must  accommodate  the  discontinuity 
in  velocity  as  the  boundary  sheet  is  crossed.  Third,  it  must  not  induce 
radial  velocities,  and  hence  radial  distortions,  at  the  inner  sheet 
surfaces.  Its  remaining  function  is  to  preserve,  in  conjunction  with 
the  flow  fields  associated  with  the  inner  sheets,  the  apparent  rigid 
axial  motion  of  the  inner  sheets  and  to  cancel  the  sum  of  vorticity 
of  the  inner  sheets.  The  first  and  second  requirements,  along  with 
the  rigid  axial  motion  constraint  and  the  constraint  on  net  vorticity, 
are  satisfied  by  the  strength  distribution  and  the  filament  geometry 
both  as  yet  unknown.  The  third  condition  may  be  automatically  satisfied 
by  a  strength  distribution  and  geometry  that  are  symmetrical  about  the 
lines  of  intersection.  This  is  simply  a  restatement  of  the  helical 
symmetry  requirements  on  the  vortex  systems. 

Having  to  solve  for  both  the  strength  distribution  and  the 
filament  geometry  of  the  boundary  vortex  sheet  presents  considerable 
difficulty  and  implies  the  need  for  some  kind  of  iterative  procedure 
for  locating  the  compatible  strengths  and  positions  of  the  boundary 
sheet  filaments.  This  difficulty  may  be  eliminated  through  a 
consideration  of  the  implications  of  the.  helical  symmetry  requirements. 


•  'V-* 
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It  has  been  established  that  the  boundary  sheet  strength  and  filament 
pitch  angle  are  constant  along  an  intersection  line.  At  an  arbitrary 
distance  Az  below  this  line,  the  strength  and  filament  pitch  angle 
will  be  different.  However,  at  any  other  angular  position  at  the 
same  Az  below  the  intersection  line  the  pitch  angle  and  strength 
must  differ  from  the  values  at  the  intersection  line  by  the  same 
amounts  due  to  the  helical  symmetry  requirement.  That  is,  at  a 
z-position  between  intersection  lines  it  is  possible  to  change  the 
sheet  strength  and  filament  pitch  angle  only  through  the  addition  of 
an  infinitesimal  strength  vortex  filament  at  the  pitch  angle  cp^. 

Thus,  it  is  possible  to  replace  the  boundary  sheet  by  two  simpler 
systems  whose  combined  effects  satisfy  all  of  the  conditions  mentioned 
above . 

The  first  of  these  sheets  is  a  uniform  sheet  of  helical  vortex 
filaments  having  constant  density  v(Cg)  a^d  constant  pitch  angle  cp^ 
as  shown  in  Fig.  6a.  This  sheet  satisfies  the  required  compatibility 
conditions  as  previously  discussed.  The  second  is  a  cylindrical  sheet 
of  helical  filaments  of  varying  and  unknown  strength  but  having  a 
constant  and  known  pitch  angle  cpR  as  shown  in  Fig.  6b.  This  sheet 
must  have  zero  strength  along  the  lines  of  intersection  in  order  to 
preserve  the  compatibility  established  by  the  first  sheet.  It  must 
have  a  symmetrical  strength  distribution  about  these  lines  and  the 
lines  midway  between  adjacent  intersection  lines  and  must  cancel  the 
radial  velocities  at  the  cylindrical  boundary  associated  with  the 
inner  sheets.  Superposition  of  the  two  sheets,  as  shown  in  Fig.  6c, 
must  maintain  the  apparent  rigid  axial  motion  of  the  inner  sheets  in 


(a)  Boundary  Sheet  of  Uniform  Strength. 


(b)  Boundary  Sheet  of  Varying  Strength. 


(c)  Result  of  Superposition. 


Fig.  6.  Concept  of  Ultimate  Wake  Vortex  System. 
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conjunction  with  the  flow  fields  associated  with  those  sheets,  and 
must  satisfy  the  net  vorticity  constraint.  These  are  the  conditions 
which  will  be  placed  on  the  solution. 

It  should  be  noted  that  to  an  observer  fixed  on  an  inner  helical 
sheet  the  boundary  vortex  sheet  (consisting  of  the  superposition  of  the 
uniform  boundary  sheet  and  the  non-uniform  boundary  sheet)  appears  in 
terms  of  the  local  strength  distribution  and. geometry  to  be  fixed 
relative  to  the  inner  sheet.  Although  the  boundary  sheet  is  actually 
slipping  forward  relative  to  the  inner  sheets,  the  observer  sees  the 
same  local  vortex  sheet  strength  at  a  given  point  on  the  boundary  at 
any  instant  while  the  identity  of  the  filaments  at  the  point  is 
constantly  changing. 

Now  consider  the  geometry  and  motion  of  the  uniform  boundary 
sheet  relative  to  the  last  outboard  filament  of  an  adjacent  inner 
helical  vortex  sheet.  The  velocity  diagram  of  Pig.  7  illustrates 
the  relationship  between  the  velocities  associated  with  the  two 
systems  according  to  the  compatibility  condition  expressed  in 
equations  (l)  and  (2).  From  Pig.  7 


V  +  w 

1  ■  ta“  ‘Pr  -  SS~ 

V  +  w_ 

\  ‘ tan  »B  =  -TST- 


(3) 


where 


* 


sin 


'Pr 


+  w  cos 


V- 


sec 


<PB 


(4) 


L 


Since  cpg  =  tan"1  X.g,  the  value  of  y(Cb)  is  known  according  to 

y(Cb)  =  *  cos  sec  (<PR  "  <Pg)  •  (7) 

The  ultimate  wake  vortex  system  of  the  heavily  loaded  ducted  fan 
is  thus  defined  in  terms  of  the  inner  helical  sheets  of  Fig.  5  (one  for 


'^^4 asasgg 


-  '.«« 
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each  blade)  having  unknown  strength  but  known  geometry,  a  vortex  sheet 
of  uniform  strength  v(£g)  8111(1  constant  filament  pitch  angle  cp^  lying 
on  the  cylindrical  wake  boundary  as  shown  in  Pig.  6a,  and  a  sheet  of 
varying  strength  but  constant  filament  pitch  angle  cp  also  lying  on 
the  cylindrical  wake  boundary  as  shown  in  Fig.  6b.  For  a  given  blade 
number  b,  pitch  and  loading  parameter  w,  the  motion  and  geometry 
of  the  system  are  determined  along  with  the  strength  of  the  uniform 
boundary  sheet.  The  solution  of  this  model  for  the  unknown  strength 
distributions  of  the  inner  sheets  and  the  non-uniform  boundary  sheet 
proceeds  directly  through  the  application  of  the  Biot-Savart  equation. 


P 


CHAPTER  III 


ANALYSIS  AND  SOLUTION  OF  THE  WAKE  MODEL 


From  the  discussion  of  the  preceding  chapter,  the  geometry  and 
motion  of  the  ultimate  wake  vortex  system  of  infinite  length  are 
known.  The  problem  is  now  to  find  the  distribution  of  vorticlty  in 
the  wake  which  will  satisfy  the  velocity  boundary  conditions  implied 
by  this  geometry  and  motion.  In  this  analysis,  the  Biot-Savart  equation 
supplies  the  required  relationship  between  the  geometry,  motion,  and 
vortex  sheet  strengths.  For  an  elemental  length  of  an  arbitrary  vortex 
filament 


The  integral  relations  for  the  velocity  components  in  Cartesian 

ip 

coordinates  are  given  by  Lamb  and  are  repeated  here  for  a  single 
finite  strength  vortex  filament. 
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A  vortex  sheet  of  the  wake  system  is  divided  into  a  number  of 
equal  width  strips.  These  strips  are  replaced  by  vortex  filaments  of 
finite  but  unknown  strengths  lying  on  the  centerlines  of  the  strips. 

The  strength  of  the  filament  must  be  equal  to  the  integral  of  the  sheet 
strength  across  the  strip  width.  For  purposes  of  numerical  calculation, 
an  adequate  representation  of  the  vortex  sheet  is  achieved  by  placing 
a  filament  along  the  centerline  of  its  corresponding  strip,  provided 
the  strip  width  is  sufficiently  small  compared  to  a  characteristic 
sheet  width. 

The  integrals  for  the  velocity  components  associated  with  the 
finite  strength  helical  vortex  filaments  introduced  in  this  manner 
are  more  conveniently  expressed  in  polar  coordinates  so  that  the 
following  transformation  is  made.  (See  Fig.  2.) 

x 7  =  r7  cos  Y7 
x  =  r  cos  Y 

y#  =  r7  sin  Y7 
y  =  r  sin  Y 

z7  =  z7  +  r'  Y7  tan  cp7 
o  T 

z  =  z 

s7  =  r7  Y7  sec  cp7 

Employing  the  transformation  yields 


filament  become 
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Aur  J  [r/2  tan  cp7  sin  (Y'  -  Y) 

—00 

+  r7  [z  -  z7  -  r7Y7  tan  cp7]  cos  (Y7  -  Y)]  ~ 

Aue  =^tJ  [{r/  taxi<pf  Cr  -  r/  cos  (*7  -  y)3 
*  -00 

+  r7  [z  -  z7  -  r7Y7  tancp7]  sin  (Y7  -  Y)|  cos  cp 
+  |r7t"  -  rr7  cos  (Y7  -  Y)j-sin  cpj 

00 

Au  =  ^  J  jjjr/2  -  rr7  cos  (Y7  -  Y)|  cos  cp 

®  —00 

-  {r7  tancp7  [r  -  r7  cos  (Y7  -  Y)] 

+  r7  [z  -  z^  -  r7Y7  tancp7]  sin  (Y7  -  Y)j-  sin  cpj 
\ 

where  the  limits  on  the  integrals  refer  to  the  infinite  extent  of  the 
ultimate  wake. 

Now,  non-dimensionalizing  these  equations  with  w  and  R,  and 
employing  the  helical  relation 

r  tan  cp  =  r 7  tan  cp 7  =  R  tan  cp^  , 
the  elemental  velocities  associated  with  one  filament  become 


-wnw^w M  i 


^U  |»«  [- 

TT  =  lr'  tan*R  sin  C"  -y) 


+  r/[z  -  z^  -  Y'  tan  (p^]  cos  (V '  -  Y)J  “■ 


^  =  ^J>  ^r(§’ 


—  / 

+  - - 2cos 

r 


(y7  -  ¥)) 


+  [z  -  z7  -  Y7  tan  cpR]  Sin  (Y7  -  Y)j  ~ 


2  r,  ?' 


-  tan  cpR  U  "  f“  cos  (Y7  -  Y)] 

-  J-  sin  (Y7  -  Y)  tan  cpR  [z  -  z7  -  Y7  tan  cpR]] 


where 


P2  =  r2  +  r/2  -  2  rr7  cos  (Y7  -Y)  +  [z-z^-Y7  tan  cpR]J 
cp  =  tan"1  [i  tan  cp  J 


and  r,  rj  z,  z7  are  non-dimensionalized  by  R. 

The  boundary  conditions  may  be  written  by  summing  the  contribu¬ 


tions  of  every  filament  of  the  system.  They  are, 


on  the  inner  sheet: 


and  on  the  cylindrical  boundary: 


(12) 

(13) 

(14) 


The  uniform  boundary  sheet  induces  no  radial  velocities  anywhere 
in  the  ultimate  wake  and  the  non-uniform  boundary  sheet  induces  no 
radial  velocities  on  an  inner  sheet  due  to  its  symmetry  above  and  below 
a  line  of  intersection.  Examination  of  the  integrand  of  equation  (9) 
reveals  that  no  radial  velocities  are  induced  at  an  inner  helical  sheet 
by  the  evenly  spaced  inner  sheets  themselves.  Thus,  equation  (13)  is 
satisfied  identically.  The  condition  on  the  cylindrical  boundary, 
equation  (l4)  involves  only  the  inner  sheets  and  the  non-uniform 
boundary  sheet.  The  condition  of  equation  (12 )  involves  all  three 
systems.  The  remaining  constraint  requires  that  the  sum  of  the 
strengths  of  all  of  the  vortex  filaments  comprising  the  wake  be  zero. 

The  boundary  conditions  may  now  be  written  on  the  inner  sheet  (from 
equation  (12)) 

^(^)^[(^X-(^^(^)-os,  (15) 


on  the  cylindrical  boundary  (from  equation  (l4)) 


30 


3 


(16) 


and 


\  (ik) +  \  Gsk) + 


Rw/ 


=  0 


(17) 


£ 


where  )  refers  to  the  inner  sheets,  )  refers  to  the  -uniform  boundary 


sheet,  )  refers  to  the  non-uniform  boundary  sheet  and 


A.  =  cos  (cp  -  cp  )  , 

B 


A  =  sin  (cp  -  cp  ) 
B 


tan  co  =  i  tan  and  tan  ®  =  —  tan  cp*, 
a  jg  r  *  j5  2*  ~x\ 


For  a  fixed  choice  of  \  and  w,  the  contributions  of  the  uniform 
boundary  sheet  to  equations  (15 )  and  (17)  are  fixed  in  terms  of  the 
sheet  strength  v(CB)  and  the  pitch  angle  cpg.  Thus  the  system  of 
equations  can  be  written  more  conveniently  as 


m 
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and 

^  (A)  +I3  (A)  =  "  \  Gife)  (20) 

Equation  (l8)  can  be  evaluated  at  a  number  of  control  points 
on  an  inner  sheet.  These  points  are  placed  between  the  filaments 
comprising  the  sheet  and  are  equal  in  number  to  the  number  of  filaments 
on  the  sheet  (the  final  point  being  placed  at  r  =  R).  Equation  (19) 
can  be  evaluated  at  control  points  on  the  cylindrical  boundary  between 
an  intersection  line  and  the  point  midway  between  intersection  lines. 

The  points  are  again  placed  between  filaments  with  one  less  point  than 
filaments.  Equation  (20 )  includes  the  filaments  on  one  inner  sheet, 
the  filaments  on  the  non-uniform  boundary  sheet  lying  between  an 
intersection  line  and  the  adjacent  intersection  line,  and  those 
filaments  of  the  uniform  boundary  sheet  passing  through  a  line 
connecting  two  adjacent  intersection  lines.  These  filaments  comprise 
a  characteristic  portion  of  the  vortex  wake,  although  the  velocities 
must  be  calculated  using  all  of  the  filaments  of  the  wake  system. 

Now,  the  integrals  of  equations  (9)  and  (ll)  may  be  defined 
respectively  as 


and  the  system  of  equations  may  be  written  as  follows: 
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for  the  control  points  on  an  inner  sheet 


\  h  (A)  0K  -  GW/~ -  (=3) 


for  the  control  points  on  the  cyli.  "ical  boundary 


\  Ir  (&)  +  ^  Zr  (A)  =  0 


(24) 


and 


\  Gflfc)  +  ^  (&)  =  -  \  Gmfe) 


(25) 


where 


and 


c 


is  over  the  filaments  of  one  inner  sheet, 

is  over  the  filaments  of  a  characteristic 
portion  of  the  non-uniform  boundary  sheet. 


is  over  the  filaments  of  a  characteristic 
w  portion  of  the  uniform  boundary  sheet. 

With  the  equations  in  this  form,  the  influence  of  w  is  confined 
to  the  right  hand  sides  of  the  equations  through  the  velocity  field  of 
the  uniform  boundary  sheet,  and  the  coefficients  Ir  and  1^  depend  only 
on  the  choice  for  A.  (the  pitch  of  the  inner  sheets)  and  b  (the  number 
of  blades). 


33 


Owing  to  the  simplicity  of  the  uniform  boundary  sheet,  the 
velocity  field  inside  the  wake  associated  with  this  sheet  has  been 
evaluated  explicity  in  terms  of  \  and  w  in  Appendix  II.  The  result 
is 


u 


jz 

w 


y(Cb) 

T- C0E  ®B 


W  W 

u2 

Using  equation  (7),  ~  can  be  written  as 


-£■  =  cos  cpR  sec  (cpR  -  cpB)  cos  cpB 
or 

u  t 

z  _  1 

w  1  +  UB  * 


Then,  the  right  hand  side  of  equation  (23)  becomes 


Further,  the  right  hand  side  of  equation  (25)  becomes 


/  v  \  y(C-d) 

(A)  =  -  TSr^  (2rTR  XB  cos  Vb) 


Since  2nR  \  /b  is  the  length  of  the  line  defining  the  character:. Stic 
portion  of  the  uniform  boundary  sheet  and  y(C_)  is  the  sheet  strength 

D 


per  unit  length.  Thus 


3b 


(27) 


Some  rearrangement  of  the  terms  of  equations  (26)  and  (27)  yields 


1  -  i/(i+xxB)  = 


(28) 


3  -  K^f)] 


(29) 


From  equation  (5) 


\  Sin  % 


W  COS  CD. 


R 


=  tan  (cpR  -  cD-q) 


B 


or 


1/  X  “  XBX 
"  X\1  +  XXR/ 


So  that,  defining 


(30) 


equations  (28)  and  (29)  may  be  rewritten  as 


1  +  X 


1  -  1/(1  +  XXB)  =  G 


H?'  Ipwws, 


£/i*  i 
|’-i 

I! 

f-'i 

K* 


II 


■rrRw/ 


.  -  £.  f-i— ) 

2b  L  .  ,27 


1  +  X 


Thus,  the  system  of  equations  can  be  written 

at  the  inner  sheet  control  points 


h  +  4  =  G 1  —  2 


at  the  cylindrical  boundary  control  points 


\  Ir  Gmfc)  +  ^  Jr  UiW  =  ° 


and 


\  (sfe)  (A) =  • 0 h (7775) 
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■1 


r 

r 


Since  the  right  hand  side  of  every  equation  is  multiplied  by  G  (or  is 
zero)  a  new  vortex  filament  strength  can  be  defined  as 

Y  =  §  (A) 

and  the  system  written 

at  the  inner  sheet  control  points 


1  +  \ 


(31) 


l  ‘0&&**i** 

•  A- 


r . 


Jjt-yi'oiiXt'M  Ui  .« 


at  the  cylindrical  boundary  control  points 


f 

V 

f 


I 

! 

I 


i 

5; 

£' 

i 

I 

.4 


{ 

l 

i- 
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<32> 

and 

\ v_+^  v'=-^(rt^)  (33) 

c  c 

The  system  of  equations  in  this  form  does  not  contain  the  w  parameter, 
so  that  a  solution  may  be  obtained  which  may  be  scaled  directly  for 
any  value  of  w.  That  is,  the  equations  are  solved  for  G  =  1  (w  =  0) 
and  the  wake  vorticity  distribution  for  any  value  of  w  is  obtained  by 
multiplying  the  result  by  the  appropriate  value  of  G. 

The  blade  bound  vortex  strength  is  found  at  any  radial  station 
by  summing  the  filaments  of  an  inner  sheet  lying  inboard  of  the  radial 
station  in  the  ultimate  wake.  The  remaining  elements  of  the  solution, 
the  strengths  of  the  filaments  of  the  non-uniform  boundary  sheet,  are 
required  in  the  calculation  of  the  power  and  thrust.  The  evaluation 
of  the  coefficients  1^  and  1^  of  the  unknown  filament  strengths,  the 
positioning  of  the  filaments  and  control  points  on  the  sheets,  and  the 
simultaneous  solution  of  the  system  of  linear  equations  for  the  unknown 
strengths  are  all  considered  in  some  detail  in  Chapter  V. 


CHAPTER  IV 


CALCULATION  OF  THRUST,  POWER,  AND  INDUCED  EFFICIENCY 

The  solution  for  the  distribution  of  vorticity  in  the  ultimate 
wake  allows  a  detailed  calculation  of  the  induced  velocities  within 
the  wake.  A  knowledge  of  these  velocities,  and  hence  of  the  momentum 
and  energy  in  the  wake,  leads  to  a  straightforward  calculation  of  the 
thrust,  power  and  efficiency. 


Thrust 

Following  the  analysis  of  Theodorsen11  consider  a  control  volume 
enclosing  the  ducted  fan  and  its  ultimate  wake  as  shown  schematically 
in  Fig.  8.  Using  the  momentum  theorem  the  thrust  of  the  ducted  fan  may 
be  found  by  considering  the  average  pressure  forces  acting  on  the  control 
surface  and  the  average  flux  of  momentum  through  the  surface.  These 
averages  are  taken  over  a  time  At  =  2rr/tG  and  the  integration  is  with 
respect  to  time,  dt  =  dz ,/(V  +  w).  Thus,  from  Fig.  8 

T  +  i  J  p  at  as  =  i  J  P  (v„  +  uj  at  <v.  +  s)  as 

SfjSj  S2 

-  ~  f  p  V  dt  V  dS 
At  1  00  00 

S1 


which  can  be  written  as 
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P 


T  +  g(v>  «)  J <*»  -  dz  dS  -  gTfljfe-5)  J <T»  “z  +  \>  42  ds  (34) 

In  order  to  integrate  equation  (3^)  the  pressure  term  must  first  be 
eliminated  by  employing  the  equation  of  motion  for  an  unsteady, 
incompressible,  potential  flow. 

H  +  Etlv2  =  f(t) 

Since  there  can  be  no  induced  velocities  at  an  infinite  radial  distance 
from  the  wake  axis 

||  +  p/p  v  V2^  =  pQ/p  (35) 

The  unsteady  term  may  be  eliminated  by  considering  the  potential 
field  in  a  steady  coordinate  system  such  that 

$(r,  Y,  z)  =  §(r,  Yq  +  Ort,  z) 


Then 


3Y 

at  "  dY  3t  "  'V 


nr 


=  X  vy 


OR 


Now  uy  may  be  written  as 


=  u-  sin  rp 
so 


cos  cp  -  u  sin 


CP 


and 


where 


(p.  -  p)  *  |  pCv2  -  vj;  )  -  p(v,  +  w)  (uz  -  ) 


V2  -  (V.  +  u. )2  +  u2  +  uj 


*  z'  r 


^  =(V„+u„  )S+4 


Substitution  of  equation  (37)  into  equation  (3*0  yields 


! - rf  fv  u  +u2-i([v  +  u  ]‘ 

+  w)  J  L  00  Z  z2\Loo  zJ 

Volume 
Inside 
Wake 


r  „  .  -.2  2  2  2  ^ 

*\**i-\) 

+  (V^  +  w)  (u^  -  u  )J  r  dr  dz  dH1 


The  limits  of  integration  are  taken  over  a  characteristic  volume  of 


the  wake  such  that 


0  s  z  s  2nRX/b 


OS  rSR 


OS  f  s  2rr 


Then  defining  non-dimensional  lengths  as 


|OS*fV*“**  " 


k 


r 


s, 

Si 

i: 


imp; — rvr^7T^rl 
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z  =  z/(2nRX/b)  ,  OS?s  1 


r  =  r/R  ,  Oi  rsl 


the  thrust  can  bo  written  as 


T  = 


a rarhj  f  J1  f  [▼.  «,  +  \  - 1  [<▼„  +  \Y 


0  0  0 


-  (V,  +U^  f  +  u2  +  -  u^  ]  +  (V,  +  v)(uz  -  u  )]r  dr  dz  dY/2n 


Now,  non-dimensionalizing  the  velocities  according  to  u  =  u/w  and 
w  =  w/QR  yields 

2  fnT>\2  f1  f1  I*2*  fl  V°°  -  -2  1  -  \2 

t  =  2«  (m)  J  J  J  Ui5^+\-2«-r-tV 

0  0  0 


,V™  .  ,2  J  .2  -  , 

-  (-S-  +  \  >  +  ur  +  -  \  ] 


Vqr 

( — = —  +  l)  (uz  -  u  )Jr  dr  dz  dY/2n 


Noting  that  (V ^QR)  =  (X  -  w)  and  defining  a  thrust  coefficient  as 

C-,  =  T/[p(nR)2TrR2] 


(33) 


the  resul t  may  be  reduced  to 
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Employing  the  definition  of  Theodorsen  for  the  non-dimensional  blade 
bound  vortex  strength  as 


K(x) 


iXk l 

2ttI<vX 


(4o) 


yields 


J  uz  dz  =  K(x)  +  ug 


so  that 


J  J  J  Sz  r  dr  dz  dY/2 n  =  J  J  [K(x)  +  u-  ]  x  dx  dY/2rr 


0  0  0 


0  0 


=  \  [2  J  K(x)  x  dx  +  u„  ]  . 


But  the  integral  term  is  defined11  as  X,  the  propeller  mass  coefficient, 


,1 

X  =  2  I  K(x)  x  dx  . 


(4l) 


Thus 


J  J  J  uz  r  dr  dz  dY/2rT  =  \  +  \  ]  (42) 

0  0  0  0 

If  r  is  taken  as  ?  =  x  =  1  then  u  becomes  u  so  that  the  integration 

Z  \ 

of  the  last  term  of  equation  (39)  becomes 
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*1  r»2n 


J  J  I  ua  r  dr  <35  di/2rr  =  ^  [K(l)  +  u^  ] 


0  C  0 


The  remaining  terms  of  the  integral  of  equation  (39) >  which 
must  be  evaluated  numerically,  may  be  defined  as 


ftl  i«X  *2rr 


pX  pX  nCA  I  p 

ez  “  J  J  J  5Z  ?  dr  d 2  d!t/2n 


0  0  0 


itX  rtl  a2tt 


n*4  pxll  p 

J  uj  r  dr  dz  dHi/2n 


0  0  0 


1*1  pi  p3rr  p 

sY  = J  J  J  Uy  r  dr  dz  d¥/2 rr 


0  0  0 


"J  J  I  [v  +  i£>«a*a,/2" 


0  0  0 


R  R 


and  e  =  e_  +  e,,  +  ew.  Employing  these  definitions  and  equations  (42) 

Z  x*  i 

and  (43),  equation  (39)  becomes 


Ct  =  w2[|x-K(X)*(|-1K  +  2  cz  -  e  +  cR  J 


It  has  been  established  that  the  flow  field  of  the  uniform 
boundary  sheet  is  known  and  that  the  vortex  sheet  strength  distributions, 
and  hence  the  flow  fields,  of  the  sheets  of  varying  strength  need  only 
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be  solved  for  the  lightly  loaded  condition.  It  is  possible,  then,  to 
calculate  the  thrust  coefficient  in  terms  of  the  w  =  0  wake  solution 
and  the  scale  factor  G.  Referring  to  equation  (39)  the  integral  terms 
can  be  modified  as  follows.  The  velocities  are  separated  into  those 
associated  with  the  inner  helical  and  non-uniform  boundary  sheets 
(variable  strength,  subscript  vs_)  and  those  associated  with  the  uniform 
boundary  sheet.  The  vs-  velocities  can  be  scaled  according  to 


=  G  u 


vs 


Further  K(x)  =  G  K(x)-  =  G  K  (x)  and  K  --  G  K  .  From 

w=o  o  o 

equation  (28)  the  velocity  associated  with  the  uniform  boundary  sheet 
GX2 

is  simply  u  =  1 - s  for  all  r.  For  simplicity  in  the  following 

Z  1  +  x2 

p  pl  pi  p2rr 

development  Jd(vol)  will  be  taken  to  represent  I  I  r  dr  dz  d¥ / 2rr . 

0  0  0 

The  integrals  in  equation  (39)  can  then  be  written  as 


Ju2  d(vol)  =  j[G  u  +  (1  -  -^-^s)]2  d(vol) 
J  z  J  Zvs  1  +  X2 


2 

=  G2  Ju2  d(vol)  +  2G  (1  -  -S-X  -g)  Juz  d(vol) 


1.  “  W  M 

+  X  vs 


+  |d-^)2  , 

1  +  X 


d(vol)  =  G2  Ju^  d(vol) 


9 


Ju^  d(vol)  =  G2  Ju2  d(vol)  , 


2 

Ju2  d(vol)  =  G2  Juz  d(vol)  +  2G  (1  -  — ~-g)  Juz  d(vol) 


1  +  X  vsr 


+  |(1„^)2  , 
1  +  \ 


Ju^  d(vol)  =  G2  J Uy  d(vol)  . 


From  equation  (42) 


Ju  a(voi)  -  |  C*„  -  -7-5I 

VS  1  +  \ 


and  from  equation  (43) 


Ju  d(vol)  =  |  [K  (1)  -  -i-j] 

VS^  1  +  \ 


Substituting  these  results  and  the  results  of  equations  (42)  and  (43) 
into  equation  (39)  yields 


CT  =  2*2  [II  (G  Ko  +  1  "  G) 


+  i  Tg2  fn2  d(vol)  +  G(1 
a  “  j  z  ' 

vs 


S  hi  '  /V  1  X 

,  ,  ,  2/  ^o  •  .  ,2' 


1  +  \ 


1  +  \ 


!tr , 


-  »•*» 


where,  for  a  given  b  and  \ ,  Kq,  Kq(i)  and  eQ  are  calculated  from  the 
w  =  0  solution.  Since  G  is  calculated  algebraically  for  a  given  \ 
and  w,  equation  (46)  provides  an  entire  family  of  values  for  whereas 
equation  (44)  provides  only  a  single  value. 


Power 

The  ideal  power  required  by  the  heavily  loaded  ducted  fan  may 
be  determined  through  a  consideration  of  the  induced  energy  loss  in 
the  ultimate  wake.  Following  Theodorsen^  the  energy  loss,  E,  is 
given  by  the  methods  of  classical  mechanics  as 


E  =  QQ  -  TV  = 
00 


Vsy  P  V,  v2  +  (|  p  V2  +  p  -  pjuz]  d  (vol)  (47) 


where  Q,  is  the  torque  of  the  fan  and  v  is  the  magnitude  of  the  total 
2  2  2  2 

induced  Velocity,  v  =  ug  +  u^  +  u^.  Substituting  from  equation  (37) 
for  (p  -  p  ),  the  energy  loss  becomes 


TVT  J  \  +  <5  VR  -wu  )]d(TOl) 


Then,  defining  a  non-dimensional  energy  loss  as 


e  =  E/[p  (OR)3  ttR2] 


and  non-dimensionalizing  velocities  and  lengths  as  was  done  for  the 
thrust  coefficient,  e  becomes 


3  pi  pl  i  °  1  2 

=  2"  J  j  j  Laz  +  f  (jf  "  1)v‘"  +  2  \  '  V  d?  d2  dV2n 

0  0  0  * 
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Again,  dividing  the  induced  velocities  into  those  associated  with  the 
uniform  boundary  sheet  and  those  associated  with  the  variable  strength 
sheets,  the  expression  for  e  can  be  calculated  in  terms  of  G  and  the 
w  =  0  solution.  All  manipulations  and  integrations  are  similar  to  those 
performed  in  the  thrust  analysis.  The  result  is  given  as 


e  .«3[o  &  +  §[l  --2^-d)  (1  --£-£>  (X  -  -i-,)  (49) 

w  n  i  ^  2.  +  \  °  1  + 


1  +  \ 


-  -14  (1  -  Jlil)  <K  (1)  -  -i--) 


2  '  2 '  'o' 

1  +  \  1  +  X  ° 


1  +  X 


+  |(|-D  (1  -  +  I  u  - 

1  +  X  1  +  X 


+  G2  (|  +  1)  fu  d(vol) 
vs 


+  G2(|-1)J[u|  +  u2  ]  d(vol) 


vs  vs 


+  G2  (1  -  -^-^fCu2  +  u!  ]  d(vol) 

1  +  X  ZySR  X 


-  2- 


~a  Ja2  52  a<vol> 


1,  «  U  V  «  *-> 

+  X  vs  vs. 


+  G' 


K  w 


"2  +  u|  ]  d  (vol)J 


vs  vs  vsr 
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Of  the  remaining  integral  terms  the  first  three  are  evaluated  in  the 

calculation  of  e  for  the  thrust  coefficient.  The  last  tvo  must  be 
o 

numerically  integrated  for  the  w  =  0  case  in  a  similar  manner. 

From  equation  (47)  the  power  can  be  written  as 

P  =  Qfl  =  TV  +  E 

CO 

so  that  defining  a  power  coefficient  as 

Cp  =  P/[p(nR)3  nR2]  (50) 

yields 

Cp  =  (X  -  w)  CT  +  e  .  (51) 


The  power  requirement  for  a  constant  diameter  wake  can  also  be 
calculated  directly  by  the  Kutta-Joukowski  theorem.  An  increment  of 
torque  is  given  by  dQ  =  bp  T  rdr  where  is  given  by 

Vial  ’  V»  +  “5  sln  ®  +  uc  cos  * 

-  \  * » (i  -  -p-i>  ■ 

x  +  X 

Then 

2 

dQ  *  p  (2rrEwX)  K(x)  [Vw  +  w  (1  -  -g-X-g)3  R2  x  dx 

?L  +  K 


and 


_  /_-\2  _  . 
=  p  yin;  ttWA 


,  \y  + 

L  oo  ' 


0 


K(x)  x  dx  -  2  G*X‘ 


<•1 

j  K(x)  -g- 

0  x 


dx 


+  X 
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CHAPTER  V 

NUMERICAL  EROCEDUEES 

In  the  preceding  chapters  the  model  and  its  mathematical  solution 
have  been  outlined  in  terms  of  the  basic  procedures  and  developments 
required.  To  obtain  such  a  solution  it  is  necessary  to  evaluate  the 
velocity  contributions  of  a  filament  of  unknown  strength  at  an  arbitrary 
location  in  the  flow  field  of  that  filament.  This  evaluation  will  be 
carried  out  by  numerically  calculating  the  integrals  of  equations  (9), 
(10)  and  (ll).  Using  these  results  a  system  of  linear  equations  will 
be  developed  in  terms  of  the  unknown  filament  strengths  of  the  system 
by  equating  the  sums  of  velocities  at  control  points  in  the  ultimate 
wake  to  the  required  normal  velocities  at  these  points.  In  order  to 
specify  the  parameters  of  these  equations,  decisions  will  be  made  about 
the  manner  of  subdividing  the  vortex  sheets  into  finite  strength  fila¬ 
ments  and  the  placement  of  control  points  on  the  sheets. 

After  the  system  of  equations  is  determined  and  solved,  the 
integrals  of  equations  (45)  and  (46)  will  be  evaluated  in  order  to 
calculate  the  thrust  and  power  of  the  ducted  fan.  In  addition  the 

values  of  K  (x),  K  and  ^  will  be  evaluated, 
o  o  0 

Evaluation  of  Velocity  Components 
The  evaluation  of  the  velocity  components  associated  with  a 
single  finite  strength  helical  vortex  filament  at  an  arbitrary  location 
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numerical  integration  of  an  integrand  which  is  solely  a  function  of  the 
geometry  of  the  filament  and  the  point  at  which  the  velocity  is  to  be 
calculated.  The  expressions  for  the  components,  equations  (9),  (lo) , 
and  (ll),  are  repeated  here  for  convenience. 

Au  .00 
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-i  /i  \ 

cp  =  tan  \=  tan  , 
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and  the  prime  1  dimensions  r' ,  z' ,  Y'  refer  to  the  location  of  an 

0 

elemental  length  of  the  vortex  filament.  The  unprimed  dimensions  refer 
to  the  location  of  the  point  at  which  the  velocity  is  to  be  calculated. 

The  computations  needed  to  numerically  evaluate  these  integrals 
may  be  simplified  by  converting  their  limits  from  -  ®  s  Y'  s  ®  to 
0  s  y'  £  oo.  This  is  done  in  the  usual  manner  by  splitting  the  integral 
at  Y7  =0,  switching  the  limits  on  the  -  «s?'s  0  portion  and 
redefining  Y7  to  -  Y7  in  the  negative  range.  The  resulting  velocity 
component  relations  are  (using  the  substitution  \  =  tan  cp^) 

0 


+  r7  [z  -  z7  -  \Y']  cos  (Y'  -  Y)]  i 

P1 

-  [r7  \  sin  (Y/  +  Y) 
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+  ?/  [X  (f“  +  ¥  "  2  G0S  ^ *  + 

-  [5  -  z7  +  XY7]  sin  (Yf  -  Y)]  —]•  dY7  , 

P2 

{O'2 -2  **'«»(*'.») 

0 

-  X2  [1  -  |^  cos  (Y7  -  Y )] 

-  X  sin  (Y 7  -  Y)  [2  -  z7  -  XY7]]  ^ 

P1 

+  £r/2  -  rr7  cos  (Y7  +  Y) 

-  X2  [1  -  t-  cos  (Y 7  +  Y)] 

+  X  sin  (Y7  +  Y)  [z  -  z7  +  XY7]]  ~  }  dY7  , 

P2 

P2  =  r2  +  r/2  -  2  rr'  cos  (Y7  -  Y)  +  [z  -  z7  -  XY7]2 
P2  =  r2  +  r/2  -  2  rr'  cos  (Y'  +  Y)  +  [z  -  z'  +  XY7]2  . 
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For  convenience  in  the  ensuing  discussion  the  integrands  of  the  velocity 
components  are  defined  by  fy>  f^  and  f  such  that 


iu 

-r'Til  *r  (*'s 


f?  <T's  r>  r'«  *>*>»' 


^  fc  «'■> *'• *• *:>  *•  «' 


The  numerical  integrations  proceed  in  a  straightforward  manner; 

Au 

-—  will  be  used  as  the  illustrative  example.  The  integration  is  perform¬ 
ed  by  subdividing  each  turn  of  the  helix  (AY '  =  2rr)  into  k  subintervals 
of  included  angle  2rr/k.  f  is  calculated  at  ¥'  =0,  2rr/k,  2(2rr/k),  ... 
k(2n/k)  and  these  values  are  stored  as  they  are  calculated.  Then  the 
contributions  of  these  segments  of  the  helical  filament  axe  summed  by 
a  trapezoidal  rule  integration  such  that 


61  (if)  =  (2rr/k)  [fr(0)  +  I  (fr(0)  +  fr(2n/k)) 

+  fr(2rr/k)  +  i  (fr(2rr/k)  +  fr(2  f-))  + 


•  •  •  +  f  (Ck  -  !]•£-) 


.  1  f  n  rri.  -  2rr  \  /. i  2r»  \'l 

g  WL*  “  AJT^  +  V*TVJ 
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(Note  that  6^  is  a  coefficient  of  velocity  subject  to  multiplication 
by  the  non-dimensional  filament  strength  y/4nRw. )  61  measures  the 

contribution  of  the  first  turn  of  the  helical  filament  (and  its 
reflection  for  -  2rr  £  Y7  £  0)  to  the  radial  velocity  associated  with 
the  filament.  This  result  is  stored  and  the  contribution  of  the  turn 
defined  by  2 nSf' s  4rc  (and  its  reflection)  is  calculated  in  the  same 
manner  yielding  62  Then  6g  is  compared  to  the  sum  of  6^  and  6g. 

A  percent  change  is  defined  as  100  x  |6g/(8^  +  6g)|  .  If  this  change 
is  not  sufficiently  small  then  63  is  calculated  and  the  comparison  is 
made  again  according  to  100  x  .  The  process  is  continued  until 


Aur 

the  percent  change  in  the  coefficient  of  is  less  than  some  specified 
tolerance.  The  coefficient  is  then  set  equal  to  the  sum  of  the  contri¬ 
butions  of  all  terms  for  which  6^  was  calculated. 

Clearly,  the  accuracy  of  this  result  depends  on  the  choice  for 
k  and  the  specified  tolerance.  The  magnitudes  of  these  parameters 
depend  in  general  on  the  geometrical  values  r,  r',  cpR  and  on  the 
number  of  blades  in  the  fan.  It  would  be  possible  to  control  these 
values  very  closely  by  a  process  of  repeated  calculations  and  compari¬ 
sons  within  a  given  integration.  For  example,  6^  could  be  evaluated 
with  k  =  90 5  then  reevaluated  with  k  =  120,  again  with  k  =  1^0,  180 
and  so  on  until  two  successive  values  were  identical  to  an  arbitrarily 
specified  number  of  significant  digits.  Then  6g  would  be  handled  in 
the  same  manner.  Similarly  the  6 ^  contributions  would  be  summed  until 
the  final  contribution  failed  to  change  the  sum  in  the  chosen  number 
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of  significant  digits.  It  is  apparent  that  such  a  procedure  could  be 
enormously  time  consuming,  and  the  choice  for  the  number  of  significant 
digits  to  be  carried  must  be  made  very  conservatively.  Choosing  the 
best  alternative  at  hand,  reliance  is  made  on  the  specification  of 
minimum  requirements  over  a  broad  range  of  cases. 

The  values  of  k  and  the  tolerances  for  the  velocity  components 
were  sized  somewhat  subjectively  on  the  basis  of  their  effects  on 
solutions  of  the  wake  vorticity  distribution.  The  method  of  sizing 
a  parameter  consists  of  fixing  all  but  one  of  these  at  very  conservative 
levels  and  systematically  varying  the  remaining  parameter  over  a  wide 
range.  This  procedure  is  repeated  for  various  values  of  b  and  X  and 
the  resulting  wake  vortex  strength  distributions  are  compared.  For 
example,  k  was  varied  over  a  range  from  90  to  720.  It  was  found  that 
for  k  greater  than  180  the  solutions  for  wake  vorticity  were  essentially 
identical  for  b  and  X  taken  over  very  wide  ranges .  That  is,  the  value 
of  Kq(1)  at  b  and  X  was  changed  by  less  than  l/2  of  one  percent  by 
increasing  k  from  180  to  720.  On  the  basis  of  this  result  k  was  fixed 
at  180  for  all  subsequent  solution,  obviating  the  need  for  any  repetitive 
calculations . 

The  permissable  tolerances  for  summing  the  6  ^  contributions  to 
A%  Auf 

— *■  and  — «•  were  treated  in  a  similar  manner.  The  minimum  values  were 
w  w 

found  to  vary  with  the  choice  of  b  and  X  but  in  general  could  be  fairly 


large  due  to  a  favorable  compensation  of  the  truncation  errors.  For 

Aur 

example,  the  value  for  — *■  for  ten  turns  of  a  filament  might  be  only 

w  Au- 

95  percent  of  the  value  for  30  turns,  but  the  errors  for  the  values 


for  the  otner  filaments  would  be  nearly  the  same  so  that  solutions  for 


wake  vorticity  (or  for  K(l))  would  differ  by  perhaps  l/k  of  one  per¬ 
du*  Auf 

cent.  Thus  it  was  decided  that  the  tolerances  on  — *•  and  — »  could  be 

w  w 

set  at  two  percent  for  sill  subsequent  solutions  without  incurring 
significant  errors. 

For  the  calculation  of  -  the  truncation  tolerance  was  found 

w 

to  be  strongly  dependent  on  the  characteristic  axial  separation  of  two 
adjacent  inner  helical  sheets  as  measured  at  the  wake  boundary.  This 
characteristic  length  may  be  non-dimensionalized  by  the  wake  cir¬ 
cumference  to  yield 


2rrR  \/b 


=  1/b 


When  z  is  small  problems  arise  in  accurately  calculating  the  radial 
c 

velocities.  That  is,  when  a  calculation  point  lies  on  the  cylindrical 
boundary,  the  radial  velocity  due  to  a  filament  in  a  helical  sheet 
lying  to  one  side  of  the  point  must  be  summed  with  the  velocity 
contribution  of  a  filament  at  the  same  radial  station  of  the  helical 
sheet  lying  to  the  other  side  of  the  point.  The  radial  velocity 
contributions  of  two  such  filaments  will  be  nearly  equal  in  magnitude 
but  large  and  opposite  in  sign.  Thus  an  tcceptable  calculation  of  the 
sum  of  these  contributions  requires  accuracy  to  a  large  number  of 
significant  digits. 

The  primary  influence  of  these  calculations  on  the  wake  vorticity 
solution  occurs  in  the  sizing  of  the  non-uniform  boundary  sheet  strength. 
If  the  percent  change  of  J  6^/^  5j  is  not  held  sufficiently  small  the 
strength  distribution  of  this  sheet  may  become  quite  irregular,  changing 
sign  and  generally  degenerating.  Fortunately,  the  magnitude  of  the 
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non-uniform  boundary  sheet  strength  becomes  very  small  (as  compared  to 

K(x))  when  'z^  becomes  small  so  that  the  effect  of  inaccuracies  in  this 

sheet  strength  is  of  secondary  importance  to  the  solution  for  K(x)  and 

to  the  calculation  of  velocity  profiles  in  the  wake.  (When *z  becomes 

zero  the  model  reduces  to  the  infinite  blade  case  for  which  the  radial 

velocities  vanish  along  with  the  strength  of  the  non-uniform  boundary 

sheet.)  It  was  found  that  the  solutions  could  be  adequately  controlled 

if,  for  Is  <  l/4  the  tolerance  was  held  to  0.1  percent;  for  1/4  <  'z 
c  c 

<l/2  the  tolerance  was  l/2  percent;  and  for  l/2  <  'z  <  1  the  tolerance 
was  held  to  1.5  percent.  Higher  values  of  are  not  considered  to  be 
of  practical  interest. 


System  of  Linear  Equations 

As  may  be  seen  in  the  preceding  section,  the  velocity  component 
contribution  of  a  single  filament  at  any  point  in  its  flow  field  is  a 
linear  function  of  the  vortex  strength  of  the  filament.  If  the 
contribution  of  n  filaments  to  a  velocity  component,  say  Aur/w,  at  a 

b 

given  point  are  summed  then  the  result  is  of  the  form 


*1  Yi  +  g2  Y2  +  S3  V3  +  •  •  •  +  gn  Yn  =  V COS  'P 


where  g^  is  the  integrated  functions  of  geometry,  f^,  related  to  the 
set  of  filaments  (one  for  each  blade)  of  strength  y^  and  to  the 
calculation  point,  r,  Y,  z.  If  the  filament  strengths  (y^)  are 
unknown  and  the  geometry  of  the  filaments  is  known,  then  specifying 
a  value  for  u^  at  the  calculation  point  yields  a  linear  equation  in 
Yl»  Yg  .  .  .  y  .  If  n  such  calculation  points  are  chosen  and  the 
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velocity  component  is  specified  at  each  point,  the  result  is  a  system 
of  n  linear  equations  in  n  unknowns.  Such  a  system  of  equations  is 
used  in  the  solution  for  the  distribution  of  vorticity  in  the  ultimate 
wake.  The  inner  helical  sheets  and  the  non-uniform  boundary  sheets 
must  be  divided  into  strips  and  the  strips  replaced  by  finite  strength 
vortex  filaments.  In  order  to  arrive  at  some  criterion  for  a  minimum 
number  of  finite  strength  vortex  filaments  to  be  used  to  represent  a 
sheet,  the  blade  circulation  for  an  optimum  free  propeller  was  calculated 
using  the  methods  of  this  chapter  for  comparison  with  the  results  of 
Goldstein^  Solutions  were  generated  using  h,  6,  8,  10,  12,  and  18 
filaments  per  sheet.  The  results,  which  are  considered  in  greater 
detail  in  Appendix  III,  indicated  that  ten  filaments  constituted  an 
adequate  approximation  to  the  vortex  sheet  (i.e.,  a  maximum  value  of 
K(x)  of  about  98  percent  of  Goldstein’s  value). 

Initial  calculations  for  the  ducted  fan  wake  model  were  performed 
using  ten  filaments  to  replace  an  inner  helical  sheet  and  eight  filaments 
to  replace  a  characteristic  portion  of  the  non-uniform  boundary  sheet. 
Because  of  the  symmetry  of  vortex  strength  in  the  boundary  sheet  the 
eight  filaments  introduce  only  four  unknowns.  Control  points  were 
placed  between  adjacent  filaments  using  ten  points  on  the  inner  sheet 
at  =0  and  three  points  on  the  cylindrical  boundary.  This  arrangement 
is  illustrated  schematically  in  Fig.  10.  At  the  ten  control  points  on 
the  inner  sheet  the  contributions  to  u^/cos  cp  =  1  are  calculated  and 
summed  to  yield  ten  equations.  At  the  three  control  points  on  the 
cylindrical  boundary  the  contributions  to  u_  =  0  are  summed  to  yield 
three  equations.  A  final  (lUth)  equation  is  written  in  terms  of  the 
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•net  vorticity.  The  forms  of  the  right  hand  elements  of  the  equations 
we  modified  to  include  the  effect  of  the  uniform  boundary  sheet  as 
discussed  in  Chapter  III  so  that  the  system  of  equations  takes  the 


*1,1  Yi  +  \f2  Y2  +  •  •  •  +  ^,14  Yi4  X  /C1  +  X  ^ 


l2,l  Y1  +  A2,2  y2 


+  •  •  •  +  *2,14  Vl4  -  X2/(l  +  X2) 


*10,1  Yi  +  *io,2  v2  +  •  •  •  +  *10,14  Yi4  "  1  /(1  +  *  ) 


*11,1  Y1  +  *11,2  Y2  + 


.  +  A, 


11,14  Yl4 


*13,1  Yl  +  *13,2  Y2  +  ’  *13,14  Yl4  0 

*14,1  Y1  +  *14,2  ?2  +  ’  *  ’  +  *14,14  Yl4  =  "  (2b)/ ^  +  X  ) 


^  through  ^j0  10  BXe  sums  ^  integrated  function,  f^,  with 

the  sum  taken  over  the  corresponding  filaments  of  each  inner  helical 
sheet.  A1  ^  through  A1Q  ^  are  the  sums  of  b  integrations  of  f^  for 
the  corresponding  filaments  of  the  characteristic  portions  of  the  non- 
uniform  boundary  sheet.  Similarly,  A^  ^  through  A^  ^0  8X6  integrations 
of  fr  for  the  inner  sheets  and  A^  ^  through  A^  ^  are  the  Integrated 
coefficients  for  the  non-unifoim  boundary  sheet.  A^  ^  through  A^10 
are  1.0  and  A^  ^  through  A^  ^4  are  2.0.  y^  through  yl0  are  the 


0 


66 

unknown  strengths,  y/ (4ttRwG)  of  the  inner  sheet  filaments  and 
through  yi1+  are  the  unknown  strengths  for  the  non-uniform  boundary 
sheet  filaments. 

Thorough  investigations  of  the  effect  of  the  total  numbers  of 

filaments  and  control  point  locations  and  their  spacing  on  the  inner 

sheets  indicate  that  the  arrangement  shown  in  Fig.  10  is  satisfactory 

and  does  not  required  modification  with *zc  (i.e.,  \  and  b).  The 

satisfactory  numbers  of  filaments  and  control  points  on  the  cylindrical 

boundary  were  found  to  be  clearly  dependent  on  *z  .  That  is,  for 

'z  <  1/4  6  filaments  and  2  control  points  suffice;  for  l/4  <  % <  1/2 
c  c 

8  filaments  and  3  control  points  suffice;  and  for  1/2  <  ~  <  1,  12 

filaments  and  5  control  points  are  adequate.  The  decisions  for 

adequate  arrangements  were  made  rather  subjectively  on  the  basis  of 

comparisons  of  a  blade  number  family  of  solutions  to  the  exact  infinite 

13 

blade  solution  by  Gray  .  A  further  criterion  which  involves  the 
smoothness  of  the  vorticity  distributions  on  the  sheets  was  used  as 
a  measure  of  the  convergence  of  the  sheet  strength  distributions. 

In  Chapter  II  it  was  noted  that  the  strength  of  the  non-uniform 
boundary  sheet  must  be  zero  at  the  lines  of  intersection  and  must 
maintain  helical  symmetry  between  lines  of  intersection.  In  order  to 
meet  these  requirements  the  sheet  strength  was  defined  as  a  series  of 
trigonometric  functions  such  that 

7  =  £  a^inPV-^s] 
n=i 
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where  N  +  1  is  the  number  of  control  points  on  the  sheet.  The  series 
is  integrated  over  a  strip  width  of  the  sheet  to  yield  the  corresponding 
finite  filament  strength.  If  2  m  represents  the  strip  width  such  that 
0)  =  1/UN  then  the  filament  strength  is 


„z  +  (0 


Y  - 


Y  dz 


U) 


“  Jl  [an  (srlri) ain  "  *> cos  [222^  "  !]] 


The  strength  distribution  of  an  inner  sheet  is  also  expressed  in 

series  form.  As  ¥  approaches  1.0  it  is  necessary  that  y  become  zero 

so  that  for  points  immediately  above  and  below  the  sheet  uf  may  vanish. 

Since  the  strength  of  the  sheet  is  y  =  -  ^  this  restriction  requires 

dr 

that  ^  be  zero  at  x  =  1.0.  Thus  the  bound  vortex  strength  may  be  written 
as 


EteL  =  “j1  B 

2nRwX  ntx  n 


,  T2n  -  1  1 

sin  I — g - tt  x  I  . 


(55) 


To  calculate ' the  strength  of  a  filament  of  the  inner  sheet  the  difference 
in  bound  strength  is  taken  at  the  radial  stations  corresponding  to  the 
edges  of  the  strip  which  the  filament  replaces.  Thus 

y  =  -  2  Bn  sin  [2-g~-  tt  ®]  cos  [2ng~  1  tt  x]  . 


The  system  of  linear  equations  is  formed  in  the  same  manner  as 
before  by  summing  the  velocity  component  contributions  at  a  given  point 
due  to  all  of  the  filaments  of  unknown  strength.  However,  ‘ihe  unknowns 
now  become  the  coefficients  of  the  two  series  for  y,  That  is,  at  the 
point  (r,  Y ,  z)  n  is  fixed  and  the  contributions  of  all  filaments  for 
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B^  (or  a„)  are  summed  to  yield  the  coefficient  A.  ..  For  example,  set 
n  n  i|ii 

n  =  1  in  the  B  series  so  the  coefficient  for  B,  it 
n  x 


The  system  formed  in  this  manner  is  equivalent  to  the  previous  system 

(with  discrete  filament  strengths)  vith  the  additional  constraints  on 

the  end  values  of  the  vorticity  distributions. 

Solution  of  the  linear  system  for  either  the  filament  strengths 

or  the  series  coefficients  is  carried  out  by  the  Gauss  method  of 

ll* 

successive  eliminations  as  presented  by  Fadeeva  .  Having  solved  for 
the  wake  vorticity  the  blade  bound  vortex  strength  distribution  may 
be  calculated  directly  from  the  series  of  equation  (55).  XQ  may  be 
calculated  by  integration  of  the  series  to  yield 


and 


Ko  = 


4b  N"1  /2n  -  1  V1 


T  n=l  V- 


•n:  (*r  -  % 


K  (x)  =  ^  Nf1  B  sin 
°  X  n=1  n  L  2 


n  x 


] 


M,q  may  be  obtained  by  a  simple  strip  integration  based  on  filament 
strength  such  that 


(*i  +  x2)4 
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Thrust  and  Power  Integrations 

In  Chapter  IV  the  thrust  and  power  coefficients  were  shown  to 
contain  integrals  of  the  velocity  distribution  through  a  characteristic 
volume  of  the  wake.  These  integrals,  which  must  be  evaluated  numerically 
are 


e 

o 


0  0  0 


vs 


Sr  "  +  +  u|  Jr  dr  dz  dY/2n 

vs  Tvs  vsR  vsR 


and  the  integrals  occurring  in  the  power  calculation 


J  J  *2  r  dr  dZ  dY/2rc 


0  0  0  vs  VSR 


and 


J1  J  J  uz  K  +  Jr  dr  dz  dY/2rr 

0  0  0  vs  VSR  VSR 


When  the  strengths  of  all  of  the  vortex  filaments  comprising  the 
ultimate  wake  are  known,  the  components  of  induced  velocity  can  be 
calculated  at  any  point  in  the  wake  by  evaluating  the  integrals  for 
each  of  the  filaments  and  multiplying  the  result  by  the  filament 
strength.  The  procedures  for  integrating  and  summing  the  contributions 
of  the  filaments  are  the  same  as  those  used  in  setting  up  the  system 
of  linear  equations.  These  results  are  then  projected  to  yield  the 
components  of  induced  velocity  in  the  r  -,  z  and  Y  -  directions. 

Due  to  the  helical  symmetry  of  the  vorticity  and  velocity  distributions 
in  the  ultimate  wake,  the  volume  integrations  may  be  performed  by 


t  kw 
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obtaining  a  detailed  knowledge  of  the  flow  field  on  a  zr  -  surface. 

This  surface  is  bounded  by  the  wake  axis,  the  cylindrical  boundary, 
an  inner  helical  sheet,  and  the  radial  line  midway  between  two  adjacent 
inner  helical  sheets.  The  surface  is  divided  into  a  network  or  grid. 

At  the  intersection  points  of  the  grid  the  induced  velocity  components 
are  calculated,  their  squares  are  calculated,  and  all  of  this  informa¬ 
tion  is  formed  into  two-dimensional  arrays  of  data.  A  typical 
integration  surface  is  shown  schematically  in  Pig.  11. 

Since  the  motion  and  vorticity  of  the  wake  are  known,  numerical 
evaluation  of  the  velocity  components  at  the  edges  of  the  Zr  -  surface 
is  not  required.  Specifically,  at  ?  =  0  uj;  =  =  0  and  ^  =  ug  .  At  z  =  0, 
uz  =  1  -  GX2/(r2  +  X2),  Uy  =  -  GX  r/(r2  +  X2),  and  ur  -  -  |  . 

At  r  =  1  ur  =  0,  uz  =  1  -  GX2/(r2  +  X2)  +  2  y  cos  cpR,  Uy  =  -  GX  r/ 

(r2  +  X2)  +  2  y  sin  cpR.  (Here,  y  is  the  strength  of  the  non-uniform 
boundary  sheet.)  At  z  =  1/2  u^  =  0  but  uz  and  Uy  must  be  evaluated 
numerically.  For  the  arrangement  shown  in  Fig.  11,  all  of  the  velocity 
components  must  be  evaluated  at  the  internal  grid  points  of  which  there 
are  27,  and  the  Uy  and  uz  values  must  be  calculated  at  the  9  lower  edge 
points . 

The  elements  of  these  two-dimensional  arrays  of  data  are  then 

2  2  2 

squared  to  form  three  new  arrays  denoted  by  u  ,  u  and  uL 

i,J  i,J  *i,J 

p 

Taking  the  u  array  as  an  example,  the  integration  may  be  perforated 

Zi,d 

by  a  simple  strip  method  such  that 
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tt 

0  o 


J  u  r  dr  dz  dY/2n  =  £  E  Z 
JQ  z  i=0  j-0  4l  J 


K 


+  \ 
i,J  1+1 J 


.  _2  ^  _2  1 

+  u  +  u 

Zi,j+1  zi+l,j+lJ 


where  (I  +  l)  is  the  number  of  grid  points  in  the  r  -  direction  and 
(J  +  l)  is  the  number  of  grid  points  in  the  z  -  direction.  For  the 

p 

network  illustrated  in  Fig.  9,  I  -  10,  J  =  4  and  4l  J  =  1600. 

For  the  integrations  involving  the  (vs)  velocities,  calculation 

of  the  array  members  is  restricted  to  the  filaments  associated  with 

the  inner  helical  sheets  and  the  non-uniform  boundary  sheet.  The  edge 

values  become:  at  r  =  0  uL  »  u  =  0,  u  =  -  l/(l  +  X  );  at 

vs  vs  vs 

2 

Z  =  0  u  =  --o  -  2^ — 2  »  "y  =  -  GX  r/Cr2  +  X2)  and  u  = 

zvs  1  +  X  r  +  X  vs  vs 

-  | r  at  r  =  1  u  =  0,  u^,  =  -  GX/(l  +  X2)  +  2  y  sin  cpR, 

c  vs  vs 

u  =  2  y  cos  cpR.  At  z  =  l/2  ur  =0  and  uz  and  u^  must  be 
Zvs  rvs  vs  vs 

evaluated  numerically. 

The  accuracy  of  the  integration  technique  is  clearly  dependent 

on  the  choice  for  I  and  J  and  on  the  strip  method  employed .  As  a  check 

on  the  method  the  u  velocity  distribution  was  calculated  for  several 
z 

values  of  X  and  b  and  integrated  using  the  method  outlined  above.  From 
equation  (42)  of  Chapter  IV 

J  uz  d  vol  =  i  [K  +  u^  ] 


Using  values  of  I  =  5  and  J  =  4  the  numerical  results  agreed  with  the 
more  exact  results  of  equation  (42)  to  within  about  1  percent  in  all 
cases.  (See  Appendix  III.)  As  a  further  check  the  values  of  e^, 
and  €  were  calculated  in  the  same  manner  for  X.  =  1.356  and  for 
increasingly  large  blade  number.  Their  convergence  to  the  exact  value 
for  the  infinite  blade  case^  was  satisfactorily  observed  and  is 
considered  in  more  detail  in  Appendix  III. 
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CHAPTER  VI 

RESULTS 

Following  the  initial  checkouts  of  the  numerical  procedures, 
the  methods  for  the  generation  and  solution  of  the  system  of  equations 
for  the  wake  vorticity  distribution  of  the  lightly  loaded  ducted  fan 
were  programmed  for  the  digital  computer.  The  system  was  solved  first 
for  a  value  of  X  =  1.356  with  b  =  2  and  4.  The  results  for  the  blade 
bound  circulation,  K(x),  were  compared  to  the  electro-potential  analogy 
results  of  Gray  .  These  comparisons  are  shown  in  Fig.  12  and  Fig.  13* 
In  both  cases  the  agreement  is  considered  to  be  good. 

Solutions  were  also  generated  for  comparison  to  the  theoretical 
2 

results  of  Tachmindji  who  used  the  velocity  potential  approach.  K(x) 
is  shown  for  the  two  methods  in  Fig.  l4  and  Fig.  15,  with  X  =  1/3  and 
2/3  and  b  =  4  in  both  cases.  The  agreement  is  again  considered  to  be 
good. 

For  comparison  to  the  exact  theory'  for  a  ducted  fan  having  an 
infinite  number  of  blades,  solutions  were  obtained  for  b  =  2,  4,  6,  8, 
12,  16,  and  24  with  X  =  1.356.  The  results  for  the  blade  bound  vortex 
strength  show  excellent  convergence  to  the  exact  solution  as  seen  in 
Fig.  16.  Since  the  measure  of  convergence  may  be  seen  more  easily  in 
terms  of  the  mass  coefficient,  the  values  of  Kq  for  increasing  blade 
number  are  shown  in  Fig.  17. 


Figure  16.  Comparison  of  the  Blade  Bound  Vorticity  for 
a  Lightly  Loaded  Ducted  Fan  for  Successive 
Numbers  of  Blades  (Digital  Computer}  with 
the  Infinite  Blade  Solution  (Exact3). 


Figure  17.  Comparison  of  Mass  Coefficients  for  Increasing  Number 
of  Blades  (Digital  Computer)  with  the  Mass  Coefficient 
for  the  Infinite  Blade  Case  (Exact3). 
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In  order  to  illustrate  the  behavior  of  the  design  and  performance 
parameters  with  the  the  variation  of  the  load  parameter,  a  family  of 
solutions  is  presented  in  Figures  18  through  21.  Fig.  18  shows  the 
basic  lightly  loaded  circulation  curves,  Kq(x) ,  for  X.  =  l/2  and 
b  =  2,  3,  4,  6,  and  8.  The  variation  of  the  scale  factor  G  with  w 
is  illustrated  in  Fig.  19.  This  result  is  plotted  versus  w/X;  since 
w  ranges  from  zero  to  X,  w/X  yields  the  convenient  range  of  zero  to 
1.0  with  w/X  =  1.0  corresponding  to  the  static  thrust  condition.  The 
variations  of  Qj,  and  Cp  with  w/X  are  shown  in  Fig.  20  and  Fig.  21  for 
the  heavily  loaded  ducted  fans  corresponding  to  the  lightly  loaded  fans 
of  Fig.  18.  The  exact  results  for  the  infinite  blade  case  are  included 
in  these  figures  for  comparison  and  to  illustrate  the  convergence  of 
the  results,  with  increasing  blade  number,  to  the  exact  solution. 

In  Appendix  IV  a  collection  of  results  is  presented  in  several 
tables.  The  range  taken  for  X  was  0  <  X  ^  1,  and  for  w/X  from  zero 
to  1.0.  For  the  values  of  X  considered,  the  blade  numbers  were  taken 
as  2,  3,  4,  6,  8,  12  and  16.  As  may  be  seen  in  the  tables,  the  higher 
blade  numbers  were  not  considered  in  the  calculations  for  the  smaller 
values  of  X.  Rather,  the  blade  number  was  increased  until  the  results 
had  closely  approached  the  exact  infinite  blade  solution. 
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APPENDIX  I 

k 

GEOMETRY  AND  MOTION  OF  THE  INNER  HELICAL  SHEETS 

The  argument  as  to  the  geometry  and  motion  of  the  wake  vortex 

pattern  of  an  optimum  ducted  single-rotation  fan  is  essentially  the 

9  3 

same  as  that  presented  by  Betz  and  Theodrosen  .  Following  these 
approaches,  consider  a  non-cptimum  ducted  fan  which  is  producing  the 
required  thrust  at  the  expenditure  of  the  necessary  amount  of  power. 

At  a  distance  behind  the  first  ducted  fan  system  such  that  the  duct 
interference  velocities  are  negligible,  arrange  a  second  ducted  fan 
having  the  same  number  of  blades  and  rotational  speed  as  the  first 
fan  and  so  phased  with  the  first  fan  that  each  blade  intercepts  one 
of  the  sheet  of  discontinuity  that  are  shed  from  the  former  fan’s 
blade  trailing  edges.  The  diameter  of  the  duct  of  this  second  system 
is  set  equal  to  the  wake  diameter  so  that  it  intercepts  the  sheet  of 
discontinuity  that  is  shed  from  the  first  duct's  trailing  edge. 

Assume  that  the  second  fan  is  mounted  on  an  extension  of  the  shaft  of 
the  first  fan  and  assume  further  that  neither  the  second  fan  nor  the 
duct  contribute  to  the  motion  of  the  wake  nor  disturb  the  flow  in 
any  way.  Similarly,  place  a  third  ducted  fan,  and  so  forth,  until  a 
large  number  of  ducted  fans  are  arranged  in  tandem,  all  mounted  on  the 
shaft  of  the  first  fan,  all  having  the  required  phase  relation,  and 
none  contributing  to  the  motion  of  the  wake  nor  to  the  thrust  or 
power  required. 
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In  general,  certain  of  the  blade  elements  of  the  first  fan  will 
be  operating  at  relatively  high  efficiencies  while  other  elements  will 

be  operating  at  relatively  low  efficiencies.  This  will  be  evident  in 

the  wake,  as  will  be  seen  later,  by  the  pitch  of  the  wake  spiral  with 
the  efficiency  being  higher  where  the  pitch  is  lower  and  vice  versa. 
Suppose  now  that  on  the  second  ducted  fan  a  positive  increment  of 
thrust  is  added  to  a  blade  element  where  the  pitch  is  low  and  an 

equal  increment  of  negative  thrust  is  added  on  the  third  fan  to  an 

element  operating  in  a  region  where  the  pitch  is  high.  The  thrust  of 
the  complete  system  remains  unchanged  but  the  third  fan  adds  more  power 
to  the  shaft  acting  as  a  windmill  than  the  second  fan  requires  to  produce 
the  thrust  increment  so  that  a  net  reduction  of  the  power  required  by 
the  system  is  realized.  (Of  course  skin  friction  is  neglected  and  it 
is  assumed  that  the  thrust  increments  are  very  small  so  that  the  power 
recovery  factor  is  100$. )  The  efficiency  of  these  added  increments  may 
be  obtained  by  considering  Fig.  22.  Using  the  Rutta-Joukowski  theorem, 
the  increment  of  thrust  is 

AT  =  pAr(0r  -  u,)  . 

The  increment  of  torque  is 

AQ  =  pAr(V  +  u  )r  . 

oo  2» 

This  gives  for  the  efficiency 

~  _  pAr(0r  - 

1  "  0  AQ  '  pAT(Vm  +  ujflr 
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The  elemental  efficiency  is  thus  simply  a  function  of  the  ratio  of  the 
apparent  velocity  of  tb»  helical  vortex  sheet  element  to  the  free  stream 
velocity.  The  process  of  adding  an  increment  of  thrust  on  a  blade 
element  of  one  fan  and  removing  the  same  amount  on  the  following  fan 
with  a  net  reduction  in  power  required  is  continued  until  no  further 
reduction  is  realized.  At  this  point  the  efficiency  of  the  added 
element  of  thrust  will  be  the  same  regardless  of  the  radius  at  which 
it  is  added.  Prom  equation  (56)  this  occurs  when  w/V^  is  the  same  at 
each  blade  station  for  the  last  fan  in  the  array.  The  wake  behind  this 
last  fan  represents  the  wake  for  the  optimum  case.  Thus  the  optimum 
condition  is  obtained  when  the  ultimate  wake  vortex  pattern  appears  to 
move  as  a  rigid  body  and  the  pitch  of  the  inner  helical  '  <  :e  spiral  is 
uniform  along  the  radius.  The  problem,  of  course,  is  the  determination 
of  the  single  ducted  fan  which  will  yield  the  same  wake  configuration 
as  the  array. 
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APPENDIX  II 


VELOCITY  FIELD  OF  THE  UNIFORM  BOUNDARY  SHEET 

The  uniform  boundary  sheet  of  the  •ultimate  wake  vortex  model 

consists  of  equal  strength  helical  vortex  filaments  of  constant  and 

equal  pitch  wrapped  on  a  right  circular  cylinder  at  the  pitch  angle  cp-g. 

The  strength  of  all  filaments  is  the  same  and  is  given  by  y,  the  sheet 

strength  per  unit  length  normal  to  the  filaments. 

This  sheet  can  be  divided  into  two  sheets  to  be  superposed. 

One  sheet  consists  of  a  system  of  infinitesimal  strength  ring  vortex 

elements  whose  axes  are  the  axis  of  the  cyli .  .irical  wake  boundary 

surface.  The  strength  of  this  sheet  is  given  by  y  cos  cp^.  The  other 

sheet  is  a  system  of  straight  line  filaments  of  infinitesimal  strength 

lying  on  the  cylindrical  boundary  surface  and  parallel  to  the  wake  axis. 

Again,  the  strength  of  all  filaments  is  the  same  and  the  sheet  strength 

per  -unit  length  normal  to  the  filaments  is  given  by  y  sin  cp_. 

B 

The  velocity  fields  of  these  two  sheets  axe  calculated  separately, 
in  terms  of  a  Biot-Savart  integration  of  the  filaments  of  the  sheets, 
and  then  superposed  to  yield  the  velocity  field  of  the  entire  uniform 
boundary  sheet. 


Element  Sheet 


Frcan  considerations  of  symmetry,  the  angular  and  radial  velocities 


associated  with  the  ring  element  sheet  are  identically  zero.  The  geo¬ 
metrical  relationsnips  necessary  to  calculate  the  axial  velocity  field  of 
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this  sheet  are  shown  in  Pig.  23.  Prom  the  figure 


dyds '  -  y  cos  cpB  RdY  dz 


h2  =  R-r  cos  Y 


h2  =  (R-r  cos  Y  )2  +  z2 


-2  2  2  2 
P  =  R  +  r  -  2rR  cos  Y  +  z 


cos  3  =  h/P 


a  =  h2/h 


so  that 


■°°  p2rr^  cos  cpB 


cos  a  RdYdz 


p®  p2rr^  cos  cpB  (h/P)  (hg/h) 

~  2  J  J  Y  5tt 


-n 


,2lT~  cos  CPB  (1  -  f  cos  Y)  dYdz 
Y  24^  * 

1  (l  +  r2  -  2r  cos  Y  +  z2) 


where  r  and  z  are  nondimensionalized  by  R.  Now  look  at  the  z- 

2 

integration  and  define,  a  =  (l  +  r  -  2r  cos  Y)  so  that 


f  dz  _  22  j®  _  2  _  , 

J  I  3/2  "  7^72  L  f  a  \l/2  *  2// 

0  (a  +  z2)  a(a  +  2^)  0  +  *0  0 


z 


Figure  23.  Geometry  of  the  Ring-Element  Sheet. 
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or 


r 


dz 


0  (a  +  z  ) 


_2.3/2 


=  2/(1  +  r  -  2r  cos  Y )  . 


This  integral  may  be  found  in  reference  15. 


Then  u  becomes 
z 


I 

*; 

.1 


\  “J 


2rr  cos 


2rr 


fl  (1  -  r  cos  Y) 

(1  +  r  -  2?  cos  Y) 


dY 


_  „  cos  <Pb  r  f* _ 

^  2tt  l  J  a  +  B  cos  Y 
0 


dY 


„2tt 


-  r 


r  cos  Y  dY  1 
J  a  +  B  cos  Y j 


where  a  =  1  +  r  and  B  =  -2r„  But 


cos  Y  _  1  B  cos  Y  _  1  f,  _ a  I 

a  +  B  cos  Y  B  a  +  B  cos  Y  5  "  a  +  h  cos  YJ  * 


so  that 


cos  Cp 


U  =  Y 

z  T 


B 


2nr 


{J 


2tt 


dY  \ 
S  +  b  cosYJ 


Now 


dY 

a  +  b  cos  Y 


1  r2"  dY 

a  J  -  B  ... 

n  1  +  sr  cos  Y 

u  a 


,  which  is  of 


the  form 


r.  2tt  ^ 

..  -  ,  which  yields  from  reference  15 

%)  X  t  cl  COS  X 
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9b 
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2n 


dx 


2tt 


1  +  a  cos  x 


(1  -  a2) 
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for  a  <  1  or,  r  <  1  (inside  the  wake). 


Then 


J 


2n 


df 

a  +  13  cos  Y 


1 

a 


2tt 


(>  -  <I)2) 


I72 


and 


cos  cpB 


v  2. 


1/2 


"1  +  r 


'1  +  r 


cos  cpB  f  lf  J2\f  2n  M  ~  COS  'PB 

r  +  H1  ‘  r  ATT-^)/  =  Y  “af“  (tt  +  n) 


=  v  — 


'(1  -  r2) 


or 


uz  =  Y  cos  cpB  . 


Straight  Line  Element  Sheet 

From  considerations  of  symmetry,  the  axial  and  radial  velocities 
associated  with  the  straight  line  element  sheet  are  identically  zero. 

The  geometrical  relationships  necessary  to  calculate  the  angular  velocity 
field  of  this  sheet  are  shown  in  Fig.  2k. 

The  velocity  at  a  point  due  to  the  filament  can  be  evaluated 
directly  from  reference  16  as 

r  ~  sin  <Pb  38 

v  =  ’  or  vi  =  y  — 5?r— 


^  sin  cpB  RdY 
=  Y - ?s== - 


»»*:«>»  rtKT'  w  *"*  **!♦*•'  ■ 
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From  the  figure  ?  =  r  +  R  -  2Rr  cos  Y 


sin  a  = 


E2  -  r2  -  E2 
2rP 


Then 


duY  ”  2rrp  ( 


r2  -  R2  +  P2\  ~ 


2rrp  \  2rP 


■j  y  sin  cpg  RdY 


so  that 


~  sin 


fB  f2n  +  r2  -  1  r2"1  dY  \ 

r  l r  r  J  a  +  £  cos  Y-f 


_2 

where  a  =  1  +  r  and  5  *  -2r.  Again,  the  integral  is  of  the  form 
,2tt 


J 

0 


dx 
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„  211/2 
(1  -  a  ) 


,  a  <  1,  or  r  <  1 


so  that 


dY 
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■  ^  (rr 


■  ^ (-/a  -  *  * ^l2) 


Thus 


~  sin  "B  far  2rrl 
"t  ‘  V  “IS-  IT  -  7!  *  0 


*■*•«< .swrtnJKd  1 


Superposition 


From  the  results  for  the  two  sheets  the  velocity  field  inside 
the  wake  due  to  the  uniform  boundary  sheet  is  given  by 


uz  =  v  cos  c pB  , 


u  =  0  , 

r  ’ 


=  0  . 
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CHECK  CASES  FOR  THE  NUMERICAL  PROCEDURES 

In  developing  the  computer  programs  to  handle  the  numerical 
calculations  required  for  the  solution  of  the  wake  vorticity  distribution 
and  the  numerical  calculation  of  thrust  and  power  coefficients,  a  number 
of  check  cases  were  developed  in  order  to  obtain  accuracy  criteria  for 
the  various  approximations  required  for  a  solution.  The  first  of  these 
checks  involves  the  replacement  of  helical  vortex  sheets  by  a  number  of 
finite  strength  vortices. 

The  solution  for  the  blade  bound  vorticity  of  an  optimum  free 
propeller  is  well-established  and  has  been  solved  for  a  number  of 
cases  by  Goldstein  .  A  numerical  model  similar  to  the  one  employed 
for  the  ducted  fan  was  developed  in  order  to  check  the  accuracy  of  the 
integration  of  the  velocity  contributions  of  the  finite  strength  helical 
vortex  filaments  and  to  establish  a  criterion  for  a  minimum  number  of 
filaments  required  to  adequately  represent  a  helical  vortex  sheet. 

The  conditions  placed  on  the  solution  were  as  follows.  The  velocity 
normal  to  a  helical  sheet  must  be  proportional  to  the  local  pitch  of 
the  sheet  along  a  radial  line  on  the  sheet,  and  the  vorticity  of  a 
single  blade  trailing  sheet  must  sum  to  zero.  The  blade  trailing 
sheets  were  divided  into  J  strips,  the  strips  replaced  by  filaments 
lying  along  the  centerlines  of  the  strips.  Control  points  at  which 
the  normal  velocity  component  was  summed  and  specified  were  placed 
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at  (J-l)  points  midway  between  the  filaments  representing  the  sheet 
oriented  at  ¥  =0.  The  arrangement  for  J  =  4  is  shown  schematically 
in  Fig.  25.  The  result  is  a  set  of  3  linear  equations  in  4  unknowns. 

The  4th  equation  is  supplied  by  the  zero  net  vorticity  constraint. 

The  system  is  illustrated  as  follows. 

^i  1  Yi  +  Ai  g  Yi  +  •  •  •  +  j  Yj  =  cos  epr 


'i  + . +  aj-x,j  *  <»s  vT 

u  -Jl 

Yx  + . +  Yj=0 


Where,  for  example,  1  is  the  sum  of  the  geometrical  integrations  of 

f^  for  the  inboard  most  filament  of  the  b  helical  sheets.  The  system 

was  calculated  and  solved  for  X  =  l/2  and  b  =2  for  values  of  J  =  4,  6, 

8,  10,  12,  and  18.  The  results  are  shown  in  Fig.  26.  Based  on  these 

results,  it  was  decided  that  the  ten  filament  system  represents  an 

adequate  approximation  to  the  vortex  sheets.  Greater  numbers  of 

filaments  improve  the  solution  very  slowly  and  it  was  felt  that  the 

correspondingly  larger  computation  times  were  not  justified. 

In  order  to  provide  a  check  on  the  calculation  and  integration 

of  the  functions  of  velocity  in  the  ultimate  wake,  the  integral  for 

u  was  set  up  in  the  manner  described  in  Chapter  V  and  evaluated  for 
z 

several  cases  of  b  at  X  =  1/2.  The  result  was  compared  to  the  exact 
value  provided  by  equation  (42)  according  to 
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J  J  J  \ f  61  **  =  I  [x  +  ^  ] 


0  0  0 


The  comparisons  are  listed  and  reduced  to  a  percent  difference  in 
Table  1,  where 


percent  difference  = 


ft 

u  d  vol  - 

rt 

lu  d  vol 

V  Z  , 

exact 

m  kj 

numerical 

[*u  d  vol 

!  exact 

As  the  table  shows,  the  difference  in  all  cases  was  of  the  order  of 
1  percent  or  less  so  that  the  calculation  and  integration  methods  were 
considered  to  be  sufficiently  accurate. 

As  an  additional  check  on  the  volume  integrations,  the  values 
of  e  ,  ev  and  e  were  calculated  for  X  =  1  and  b  =  2,  4,  8,  and  12. 
These  results  are  compared  to  the  values  for  an  infinite  blade  fan. 

As  can  be  seen  in  Fig.  27  the  results  show  satisfactory  convergence 
for  increasing  blade  number. 


103 


Table  1.  Accuracy  of  the  Wake  Integrations 


b  Juz  d  vol  Juz  d  vol  Difference  $  Difference 


numerical 

exact 

2 

•  520 

.526 

.006 

l.l4 

3 

•  552 

.557 

.005 

0.90 

4 

.569 

•  572 

.003 

0.53 

8 

.589 

•  590 

.001 

0.17 

X  =  1/2,  v  =  0 
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Like  the  free  propeller  in  axial  flight,  a  single-rotation  ducted  fan  of  highest 
induced  efficiency  is  characterized  by  an  ultimate  wake  vortex  system  shed  from  the 
blade  trailing  edges  whose  apparent  motion  is  that  of  rigid  helical  surfaces.  In 
addition  and  concentric  with  this  inner  sheet  system  there  is  a  cylindrical  surface 
•  of  helical  vortex  filaments  shed  from  the  duct  training  edge .  For  zero  hub  diameter 
and  neglecting  compressibility,  viscosity  and  tip  clearance,  a  consistent  mathemati¬ 
cal  model  of  the  constant-diameter  vortex  wake  is  developed  and  the  compatibility 
relationships  to  be  satisified  are  presented.  Using  the  Biot-Savart  equation,  the 
vortex  strength  distribution  in  the  ultimate  wake  is  determined  and  then  related  to 
the  blade  bound  vortex  strength  distribution.  In  addition,  expressions  are 
developed  for  the  thrust,  power  and  induced  efficiency  which  depend  on  numerical 
integrations  of  velocity  profiles  in  the  ultimate  wake.  It  is  shown  that  the  wake 
vorticity  and  the  veloc  ty  distribution  in  the  wake  for  all  loadings  from  the  lightl; 
loaded  limit  to  the  heavily  loaded  static  thrust  condition  may  be  extracted  from  the 
lightly  loaded  result  through  the  use  of  a  simple  algebraic  scale  factor.  The 
performance  parameters,  thrust  power  and  blade  bound  vortex  strength  are  thus 
expressed  in  terms  of  a  lightly  loaded  solution  for  a  given  freestream  velocity, 
blade  tip  speed  blade  number  and  the  loading  parameter.  The  lightly  loaded  case  is 
compared  to  existing  theoretical  and  experimental  results.  Some  sample  results  for 
heavily  loaded  performance  parameters  and  design  parameters  are  presented  and  comparet 
to  the  exact  results  for  a  heavily  loaded  ducted  fan  with  an  infinite  number  of  blades 
Results  for  the  heavily  loaded  system  over  a  broad  range  of  conditions  are  presented. 


DD  ,'r..1473 


REPLACES  DD  FORM  147*.  I  JAN  44,  WHICH  IS 
OBSOLETE  FOR  ARMY  USB. 


Unclassified 


nrlty  Classification 


